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Abstract. The paper discusses recent progress in understanding statistical 
properties of eigenvalues of (weakly) non-Hermitian and non-unitary random 
matrices. The first type of ensembles is of the form J = H — iT, with H being 
a large random N X N Hermitian matrix with independent entries "deformed" 
by a certain anti- Hermitian N X N matrix if satisfying in the limit of large 
dimension A'^ the condition: TvH^ oc A'^TrF^. Here V can be either a random or 
just a fixed given Hermitian matrix. Ensembles of such a type with F > emerge 
naturally when describing quantum scattering in systems with chaotic dynamics 
and serve to describe resonance statistics. Related models are used to mimic 
complex spectra of the Dirac operator with chemical potential in the context of 
Quantum Chromodynamics. 

Ensembles of the second type, arising naturally in scattering theory of 
discrete-time systems, are formed hy N X N matrices A with complex entries 
such that A^A = I — T. For T = this coincides with the Circular Unitary 
Ensemble, and < T < I describes deviation from unitarity. Our result amounts 
to answering statistically the following old question: given the singular values of 
a matrix A describe the locus of its eigenvalues. 

We systematically show that the obtained expressions for the correlation 
functions of complex eigenvalues describe a non-trivial crossover from Wigner- 
Dyson statistics of real/unimodular eigenvalues typical for Hermitian/unitary 
matrices to Ginibre statistics in the complex plane typical for ensembles with 
strong non-Hermiticity: < Tr_ff^ >oc< TrF^ > when — > oo. Finally 
we discuss (scarce) results available on eigenvector statistics for weakly non- 
Hermitian random matrices. 



1. Introduction 

As is well-known, the statistics of highly excited bound states of closed quantum 
chaotic systems of quite different microscopic nature is universal. Namely, it turns 
out to be independent of the microscopic details when sampled on energy intervals 
large in comparison with the mean level spacing A, but smaller than the so called 
Thouless energy scale. The latter is related by the Heisenberg uncertainty principle 
to the relaxation time necessary for the classically chaotic system to reach equilibrium 
in phase space Moreover, the spectral correlation functions turn out to be 

exactly those which are provided by the theory of large random Hermitian matrices 
with independent, identically distributed Gaussian entries. The correspondence 
holds in the limit of large matrix dimension on the so-called local scale. The 
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local scale is determined by the typical separation A ~ {Xi — Xi^i) between 
neighbouring eigenvalues situated around a point X, with the brackets standing 
for the statistical averaging. Microscopic arguments supporting the use of random 
matrices for describing the universal properties of quantum chaotic systems have been 
provided recently by several groups, based both on traditional semiclassical periodic 
orbit expansions El and on advanced field-theoretical methods Ej ■ 

In parallel, it was realized that another type of random matrices - those 
with the so-called chiral structure - serve as a pattern for universal statistics of 
low- lying eigenvalues of the Euclidean Dirac operator. The latter is of crucial 
importance for describing such a fundamental phenomenon as spontaneous breakdown 
of chiral symmetry in Quantum Chromodynamics (QCD). That line of development 
is thoroughly described in [S]. The corresponding (anti-)Hcrmitian matrices are of the 
form 

where J stands for a complex matrix, with being its Hermitian conjugate. Such an 
off-diagonal block structure is characteristic for systems with chiral symmetry. One of 
the main objects of interest in QCD is the so-called partition function used to describe 
a system of quarks characterized by Np flavours and quark masses m f interacting with 
the Yang-Mills gauge fields. At the level of Random Matrix Theory the true partition 
function is replaced by the matrix integral: 

ZNFim)= dJdP Y[det{D + mf} e-^^''^^'^ (2) 
f=i 

Here instead of an integration over gauge fields in the QCD partition function we have 
the integration over matrices J It represents, in fact, the ensemble average of the 
product of characteristic polynomials of D. 

The Hermitian ensemble of block off-diagonal matrices —iD with independent, 
identically distributed Gaussian entries is known as the chiral Gaussian unitary 
ensemble (chGUE). The eigenvalues of chiral matrices appear in pairs ±Xk , k = 
1, ...,N. It is evident that the nonzero ea; = A| 7^ are the singular values of the matrix 
j, i.e. the eigenvalues of J^J. The latter are described by the so-called Laguerre 
ensemble which independently emerged in the description of quantum transport in 
mesoscopic systems jZj, see e.g.|Hj and references therein. 

All these facts make the theory of random Hermitian matrices of various kinds 
a powerful and versatile tool of research in different branches of modern theoretical 
physics PE]. 

Very recently complex eigenvalues of non-Hermitian random matrices have 
also attracted much research interest due to their relevance to several branches of 
theoretical physics. Most obvious motivation comes from the quantum description 
of open (or scattering) quantum chaotic systems EH HH IT^ whose fragments can 
escape, at a given energy, to infinity or come from infinity. For systems of this kind the 
notion of discrete energy levels loses its validity. Actually, chaotic scattering manifests 
itself in terms of a high density of poles of the scattering matrix placed irregularly 
in the complex energy plane. Each of these poles, or resonances, £k = Ek — fTfe, is 
characterized not only by energy but also by a finite width Tk defined as (twice) the 
imaginary part of the corresponding complex energy and reflecting the finite lifetime 
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of the states in the open system. NaturaUy, such a picture stimulated various groups 
to develop a statistical description of resonance parameters ^1 El El E| ■ 

Recently, the progress in numerical techniques and computational facilities made 
available high accuracy patterns of resonance poles for realistic atomic and molecular 
systems in the regime of quantum chaos, see e.g. ^| El 1201 1^ as well as for 
such model quantum chaotic systems as quantum graphs |22j and a quantum particle 
under time- and space-periodic perturbation j28| . Due to the presence of resonances 
the elements of the scattering matrix show irregular fluctuations as functions of the 
energy of incoming waves, see |24| and references therein. The main goal of the theory 
of quantum chaotic scattering is to provide an adequate statistical description of such 
a behaviour. 

The most natural framework for incorporating a random matrix description of the 
chaotic scattering (and hence for addressing statistics of resonances) is the so-called 
"Heidelberg approach" suggested in the pioneering paper and described in much 
detail in ^2 El- The starting point of that approach is a particular formulation 
of the scattering theory based on the notion of an "effective Hamiltonian" . Such 
formulation goes back to the ideas introduced in the context of Nuclear Physics 
in classical works by Kapur and Peierls, Wigner, Bloch, Feshbach and others. It 
was further developed, refined and clearly presented in the book by Mahaux and 
Weidenmiiller , see also recent discussions in |E1 El ■ Following that approach 
one introduces a division of the Hilbert space of the scattering system into two parts: 
the " interaction region" and the " channel region" . The channel region is supposed to 
describe a situation of two fragments being apart far enough to neglect any interaction 
between them. Under these conditions their motion along the collision coordinate is 
described by a superposition of incoming and outgoing plane waves with wavevectors 
depending on the internal quantum states of the fragments. We assume that at given 
energy E exactly M different quantum states of the fragments are allowed, defining 
M "scattering channels" numbered by the index a. At the same time, the second 
part of the Hilbert space is to describe the situation when fragments are close to one 
another and interact strongly. Correspondingly, any wavefunction of the system can 
be represented as a two-component vector, with first (second) component describing 
the wave function inside the interaction (channel) region. 

Using the standard methods of scattering theory exposed in detail in El El 
one can relate the two parts of the wave function to one another and finally arrive at 
the following representation of the energy-dependent scattering matrix S in terms of 
an effective non-Hermitian Hamiltonian Tie// = H — iV: 

Sab{E) - Sab - 2zn ^ [E - Ueff^^W,, 

with the Hermitian Hamiltonian H describing the closed counterpart of the open 
system (i.e. interaction region decoupled from the channel one) and the anti-Hermitian 
part iV arising due to a coupling to open scattering channels. In this expression the 
Hamiltonian H is written in some arbitrary basis of states such that Hij = {i\H\i). 
The amplitudes VFoi, a — 1,2,...,M are matrix elements coupling the internal 
motion in an "internal" state \i) to one out of M open channels a. One also has 
to choose the anti-Hermitian part to be = T^^a^io-W*^ in order to ensure the 
unitarity of the MxM scattering matrix S{E) [211 ■ The eigenvalues £k of the operator 
Tie/ f = H — iV are poles of the scattering matrix and have a physical interpretation 
as resonances: long-lived intermediate states to which discrete energy levels of the 



Random matrices close to Hermitian or unitary 



4 



closed system are transformed due to coupling to continua. In particular, one can 
show tliat|ll| 

det(E -nlff) N (w f*\ 

expressing the total scattering phaseshift </>(£') = logdetS'(£') in terms of the 
egenvalues £k- 

Applying these constructions to the quantum chaotic scattering it is natural 
to expect that universal properties of scattering systems are inherited from the 
corresponding universality of the energy levels of their closed counterparts. Thus, 
we expect the scattering characteristics (in particular, the statistics of resonances and 
phaseshifts) could be adequately described by a random matrix approximation on 
a scale comparable with the typical level spacing A. In contrast, a smooth energy 
dependence of 5'-matrix elements on a much larger energy scale must be system- 
specific. 

According to the general idea one incorporates the random matrix description 
of quantum chaotic systems by replacing the Hamiltonian iJ by a random matrix 
of appropriate symmetry. Namely, chaotic systems with preserved time-reversal 
invariance (TRI) should be described by matrices Hij which are real symmetric. 
On the other hand systems with broken TRI are modelled by complex Hermitian 
matrices from the Gaussian Unitary Ensemble (GUE) 1 . The last step of the 
Heidelberg approach is performing the ensemble averaging non-perturbatively. Usually 
it is done in the framework of the so-called supersymmetry technique introduced 
initially by Efetov [221 and adopted to the description of quantum chaotic scattering 
by Verbaarschot, Weidenmiiller and Zirnbauer 

The Heidelberg approach turns out to be a very powerful tool for extracting 
different universal properties of open chaotic systems, and many important quantities 
characterizing the chaotic scattering were successfully investigated along these lines 
p.lj, 11 . 12 .3Qt . In particular, analytical predictions of the distribution of resonance 
widths jlll ITS] were found to be in good agreement with available data obtained 
numerically for realistic models of chaotic quantum scattering \22i 1231 15T] . 

In fact, the outlined description of the wave scattering can be looked at as an 
integral part of the general theory of linear dynamic open systems in terms of the 
input-output approach. These ideas and relations were developed in system theory 
and engineering mathematics many years ago, see e.g. j77j going back to the pioneering 
works by M. Livsic [2HI- Livsic himself stressed equivalence of his mathematical 
constructions to those used by nuclear physicists, but that theory and subsequent 
developments went practically unnoticed by the quantum-chaos community. A brief 
description of the main constructions and interpretations of the linear open systems 
approach (in particular, a short derivation of the equation Eq.Q) can be found in 
|53i| and is not repeated here. Let us only mention that the unitary scattering matrix 
S{E) is known in the mathematical literature as the characteristic matrix-function of 
the non-Hermitian fundamental operator Tigf f . 

In the scattering theory leading to Eq.© the time was considered to be a 
continuous parameter. On the other hand, a very useful instrument in the analysis of 
classical Hamiltonian systems with chaotic dynamics are the so-called area-preserving 
chaotic maps 32, 33 . They appear naturally either as a mapping of the Poincare 
section onto itself, or as result of a stroboscopic description of Hamiltonians which are 
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periodic functions of time. Their quantum mechanical analogues are unitary operators 
which act on Hilbert spaces of finite large dimension N , and are often referred to as 
evolution, scattering or Floquet operators, depending on the given physical context. 
Their eigenvalues consist of N points on the unit circle (eigenphases) . Numerical 
studies of various classically chaotic systems suggest that the eigenphases conform 
statistically quite accurately the results obtained for unitary random matrices of a 
particular symmetry (Dyson circular ensembles). 

Let us now imagine that a system represented by a chaotic map (" inner world" ) is 
embedded in a larger physical system ("outer world") in such a way that it describes 
particles which can come inside the region of chaotic motion and leave it after some 
time via M open channels. Models of such type appeared, for example, in ^J. 

The general linear system theory |2Z] considers dynamical systems with discrete 
time as frequently as those with continuous time. The corresponding construction is 
discussed in |^ and we mention here only its gross features. For a closed linear system 
characterized by a wavefunction 5" the "stroboscopic" dynamics amounts to a linear 
unitary map ^'(n) — > ^1^(71 + 1), such that \l/(n + 1) = u\E'(n). The unitary evolution 
operator u describes the "closed" inner state domain decoupled both from input and 
output spaces. Then a coupling that makes the system open must convert the evolution 
operator u to a contractive operator A such that 1 — A^A>Q. It is easy to understand 
that one can always choose A — u^l — tt^ where the matrix f is a rectangular 
N X M diagonal with the entries Tij = dijTj , 1 < i < N , 1 < j < M (forM < 
N), < Tj < 1. Such a form automatically ensures 1 - A^A = rrt > and 
suggests a clear interpretation of the constituents of the model. Namely, the equation 
\E'(n + 1) = A\['(n) describes an irreversible decay of any initial state ^"(0) 7^ when 
an input signal is absent. On the other hand, assuming a nonzero input and zero 
initial state ^'(0) = one can relate the (discrete) Fourier-transforms of the input and 
output signals at a frequency uj to each other by a M x Af unitary matrix S{uj) given 

by 

S(lu) = Vl — f I'f — f ^ -iif . (5) 

e-'^^l-A 

This equation is a complete discrete-time analogue of Eq.©. In particular, using the 
idenity det ^1 — ^1^2^ = det ^1 — ^2^1^ vahd for two arbitrary matrices Ai^2 as 
well as the relation , ^ _ = —^== it is straightforward to verify that 

dot fit _ e''^) N / * _ 

det^(.) = e-- ) . = e-- ^ ^g-^, (6) 

det le~"^ - Aj k=i ^ ~ ^' 

where Zk stand for the complex eigenvalues of the matrix A which all belong to the 
interior of the unit circle: \z\ < 1. This relation is an obvious analogue of Eq.Q and 
gives another indication that z„ indeed play the role of resonances for the discrete 
time systems. 

Now assume that the motion outside the inner region is regular. According to 
the general strategy one should be able to describe generic features of quantized maps 
representing chaotic inner dynamics by choosing the matrix u to be a member of the 

Dyson circular ensemble of random matrices. Then one finds: t^t = 1 — [(^(a;)!^, with 
the bar standing for the averaging of S'(w) in Eq.Q over it. Comparing this result 
with |10| il5| we see that the M eigenvalues < < 1 of the M x M matrix r'^'f 
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play the role of the so-called transmission coefficients and describe a particular way 
the chaotic region is coupled to the outer world. 

This line of reasoning was motivated by recent papers reviewed in j23| . The 
authors considered the Floquet description of a Bloch particle in a constant force and 
periodic driving. After some approximations the evolution of the system is described 
by a mapping: c„+i = Fc„, where the unitary Floquet operator F = Stj is the product 
of a unitary "M-shift" S : Ski = Si^k-M, l,k = — oo, ...,oo and a unitary matrix U. 
The latter is effectively of the form U = diag(di, u, (^2), where di^2 are (semi)infinite 
diagonal matrices and u can be taken from the ensemble of random N x N unitary 
matrices. One can check that such a dynamics can easily be brought to the standard 
form discussed above, with all M diagonal elements of the N x M matrix r are equal 
to unity I . 

A direct inspection immediately shows that the non- vanishing eigenvalues of the 
fundamental operator A for the latter example coincide with those of a {N — M) x 
{N — M) subblock of the random unitary matrix u. Complex eigenvalues of such 
"truncations" of random unitary matrices were studied analytically by Zyczkowski 
and Sommers in and we reproduce their study in Sec. 3.1 . 

In fact, truncated unitary matrices represent only a particular case of random 
contractions A. The ensemble of general N x N random contractions A = uVl — tt^ 
describing a chaotic map with broken time-reversal symmetry can be shown to have 
the following probability measure in the matrix space: 

V{A)dA cx 5{A^A - G)dA , G = 1 - f (7) 

where dA = J^^^- dRs Aijdlm. Aij and the normalization constant is independent of G. 

It follows by averaging the matrix 5-function 5{A — ii^/l — ff t) over Dyson's circular 
unitary ensemble. The N x N matrix ff ^ = 1 — G > is natural to call the deviation 
matrix and we denote it T. It has M nonzero eigenvalues coinciding with the transmis- 
sion coefficients Ta introduced above. The particular choice T^^^^j^j) = 1, T^^yj^j) = 
corresponds to the case considered in [21] . Eq. lO describes an ensemble of subunitary 
matrices A with given eigenvalues of A''' A, i.e. with specified singular values of the 
matrix A. The question of characterizing the locus of complex eigenvalues for a matrix 
with prescribed singular values was considered in classical papers by Horn and Weyl 
|37|. In essence, we provide a statistical answer to this question. Indeed, rotating the 
matrices A by a general unitary transformation and averaging over the unitary group 
yields a natural measure on matrices A with given singular values. In section 3 we will 
derive the corresponding joint density of eigenvalues and their correlation functions. 

Although the initial interest in non-Hermitian and non-unitary deformations of 
random matrices had its origin in the theory of quantum chaotic scattering, soon 
after it got a much of stimulus from the QCD applications of chiral random matrices. 
The reason is that in order to describe Dirac eigenvalues in the presence of non-zero 
chemical potential one has to replace the (anti-)IIermitian Euclidean Dirac operator, 
see Eq.(^, with an operator of the following structure: 

J Actually, the original paper employed a slightly different but equivalent construction dealing 
with an "enlarged" internal space of the dimension N -j- M. We prefer to follow the general scheme 
because of its conceptual clarity. 
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Such matrices are no longer (anti-)Hermitian and their eigenvalues go away from the 
(imaginary) axis to the complex plane. This phenomenon clearly bears some similarity 
to the effect of opening a closed system by converting its eigenvalues to resonances. As 
the result, the fermion determinant in the partition function Q has a complex phase. 
On the level of true lattice simulations such a phase results in major computational 
difficulties. All this makes random matrix models a very useful tool for getting insights 
into gross universal features of the analytical structure of the QCD partition function. 
In turn, information on such a structure is based on the knowledge of statistics of 
complex eigenvalues of the Dirac operator with non- vanishing chemical potential which 
stimulated a lot of research of such eigenvalues ,,38^ .39_, JO , 411 l42[B5ll44| . 

Finally, it deserves to be briefly mentioned that objects related to weakly non- 
Hermitian random matrices emerge also in the context of interesting problems of 
motion and localization of a quantum particle in disordered media subject to an 
imaginary vector potential |45| or a random scalar potential |46| . The first of these 
problems is related to the motion of superconductor flux lines in a sample with 
columnar defects, the second to the motion of a self-repelling polymer chain and 
other interesting problems (see the cited references). These developments clearly show 
that, apart from being a rich and largely unexplored mathematical object, weak non- 
Hermitian deformations of random matrices enjoy direct physical applications and 
deserve to be studied in more detail. 

The goal of the present paper is to give an overview of results and methods 
for dealing with spectra and eigenvectors of various non-Hermitian and non- 
unitary deformations of classical ensembles of random matrices. For the reasons 
which should be clear from our discussion above we concentrate on the regime of 
weak non-Hermiticity recognized in |47| . It occurs when the imaginary parts of 
typical eigenvalues are comparable with the mean separation between neighbouring 
eigenvalues along the real axis. A simple analysis shows that the last condition 
is satisfled when the traces of the Hermitian matrix H and the perturbing anti- 
Hermitian matrix iT are related as: limyv-too ^ ""^^^ discuss at 

all interesting questions related to strongly non-Hcrmitian matrices, see e.g. papers 
| 48| and references therein. 

The first weakly non-Hermitian ensemble analysed in full generality j50| was 
that of non-Hermitian Gaussian deformations, with the anti-Hermitian part taken 
independently from the GUE. For this case one can develop a rigorous mathematical 
theory based on the method of orthogonal polynomials. Such a method allows one 
to study correlation properties of complex spectra to the same degree as is typical 
for earlier studied classes of random matrices 0^1 • We briefly discuss the results 
obtained in this model in subsection 2.1 (the detailed exposition of the method and 
the derivation of the results can be found in ^SU|) as well as recent progress for other 
Gaussian ensembles and closely related models. We also introduce an alternative 
way of addressing the problem which places the correlation function of spectral 
determinants as primary object. This second method turns out to be indispensable 
when one deals with non- Gaussian deformations as those emerging in the theory 
of chaotic scattering. The latter problem was successfully studied in the work by 
Fyodorov and Khoruzhenko I5T1 and we discuss in subsection 2. 2 the main ingredients 
of their solution. At the end of that subsection we indicate a close similarity between 
the results obtained for the Gaussian and for the finite-rank deformations. 

Unfortunately, it is not known at the moment how to extend the exact methods to 
the case of non-Hermitian deformations of real symmetric matrices. The last question 
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is of great practical interest for the application of the obtained results to chaotic 
scattering. Indeed, the majority of quantum chaotic systems possesses time-reversal 
invariance and hence should be adequately described by random matrices which are 
real symmetric PP. Restricting our attention to the simplest, still highly non-trivial 
statistical measure - the mean density of complex eigenvalues of the resulting complex 
symmetric effective Hamiltonian Tie/ / -the answer can be given in a closed form |52j . 
The object turns out to be accessible via Efetov's supersymmetry approach which we 
shortly discuss in subsection 2.3. 

Very recently the properties of bi-orthogonal eigenvectors of non-Hermitian 
random matrices attracted considerable interest of various researchers [52i iMj 55 , 5B| 
IS7| . In particular, an essential boost to these studies was given by the demonstrated 
relevance of the non-Hermitian effective Hamiltonian T-Lef f = H — iV to the problem 
of lasing in random media EHl 11^1 ■ Indeed, for systems of such a type the non- 
orthogonality of eigenmodes is known to play an important role [001 ■ Similar effects 
were discussed in relation to the unimolecular decay in NO2 molecules ^J. In section 
4 we discuss this issue for the regime of weak non-Hermiticity. 

Most of our own results discussed in this review were reported earlier in the form 
of short communications [HHl IHl IK^ 1571 IHTl 1^ IH^ and in an unpublished paper by 
one of the authors jOl]. Some (but by far not all) details of the underlying derivations 
are elucidated in the present text, as far as space restrictions allow. In particular, 
appendix A contains the derivation of the formulas related to various aspects of the 
Itzykson-Zuber-Harish-Chandra integration over the unitary group extensively used 
in both [HI] and [HSI- 



2. Non-Hermitian Deformations of Random Matrices 

2.1. Gaussian deformations: from Wigner-Dyson to Ginibre eigenvalue statistics. 

In the present section we concentrate on a particular case of almost-Hermitian random 
matrices with i.i.d. entries J = H + ivA where H and A are taken independently 
from the Gaussian Unitary Ensemble (GUE) of Hermitian matrices X = with 

probability density V{X) = Q^^ cxp - Nl2J'^ Tr . 

Let us now introduce a new parameter r = {1 — v^) / {1 + v^) and choose the scale 
constant Jg to be equal to (1 -I- t)/2, for the sake of convenience. The parameter r 
controls the magnitude of correlation between Jjk and Jkj'. (JjkJkj) = t/N, hence the 
degree of non-Hermiticity. This is easily seen from the probability density function 
for our ensemble of the random matrices J: 

V{J, J^) = Cjv' exp (^-j^^-^ Tr(Jjt _ r Re P)^ , (9) 

where Cjy — [n^{l — t^)/A^^]^^/^. All the Jjk have zero mean and variance 
(I Jjfep) = and only Jjk and Jkj are pairwise correlated. If r = all the Jjk are 
mutually independent and we have maximum non-Hermiticity. When r approaches 
unity, Jjk and J^j are related via Jjk ~ Jkj and we are back to an ensemble of 
Hermitian matrices. 

Our first goal is to determine the n-eigenvalue correlation functions in the 
ensemble of random matrices specified by Eq. The density of the joint distribution 
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of eigenvalues in this ensemble is given by 

-N 



pAr({Z}) = CAr(T)|A({Z})|^exp 



N 



1 -r2 ^ 



[\Z, 



Zf) 



(10) 



where cn{t) 



jYN(N + 1)/2 



{Z} = (Zi, . . . , Zjv) is the vector of eigenvalues and 



7i-«l!---Ar!(l-T2)»/2 • 

A({Z}) = YijKki^j ~ is the corresponding Vandermonde determinant. 

To derive Eq. H1U|) we integrate 'P{J,J^) from Eq. @ over the surface of all 
complex matrices whose eigenvalues are Zi,...Zn. Following Dyson (01], p. 501, 
see also jtiSI ) we decompose every complex matrix with distinct eigenvalues as 
J = U{Z + R)U\ where Z = diag{Zi, . . .Z^}, C/ is a unitary matrix, and -R is a 
strictly upper-triangular one. If we label the eigenvalues and require the first non-zero 
element in each column of U to be positive, then the decomposition is unique (it is 
frequently referred to as "Schur decomposition"). The Jacobian of the transformation 
J — » {Z,R,U} depends only on Z and is given by the squared modulus of the 
Vandermonde determinant. So, integrating out R and U is straightforward and the 
resulting expression is Eq. (^UJ- 

The correlations between eigenvalues are conventionally characterized by the n- 
eigenvalue correlation functions. To define them we subdivide the vector of eigenvalues 
{Z} = {Zi,...,Zn) into two parts {z} = (Zi,...,Z„) and {^} = (^i, ...,^fc, ...,^Ar_„) 
identifying = Zn+k- Then the correlation functions are given by: 

M{z}) - 777^ / d'{^} 'Pn{{z}, m (11) 



{N- 



where d"^ {S} = Y{^S^^ dKe dlm^/. 

The form of the distribution Eq. (|10() allows one to employ the powerful method 
of orthogonal polynomials 49 . Let Hn{z) denote the n-th Hermite polynomial, 



{±ir 

V2^ 



exp I ^ 



dt exp [ ^izt 



(12) 



The crucial observation borrowed from the paper (see also the related paper [S7]) 
is that the polynomials 




n = 0,1,2,, 
function 



Pn{Z) = 



are orthogonal in the complex plane Z = X 



(13) 

iY with the weight 



w'^{Z) = exp 



N 



(l-r2) 



\Z 



-{Z^ + z* 



i. e. J d:^Zpn{Z)p,n{Z*)w^{Z) = (5„™, where d^Z = dXdY. 

The standard machinery of the method of orthogonal polynomials 49 yields the 
functions Rn{{z}) in the form 

i?„({z}) = det [KNiZ,, ZD] \j,k=i,...,n , (14) 

where the kernel Kisi{Zi, Zj) is given by 



N-l 



Kn{Zi,Z2) = w{Zi)w{Z2) 22 Pn{Zl)Pn{Z2) 

n=0 



(15) 
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With Eqs. (|T3J|-||I3) at hand, let us first examine the regime of strong non- 
Hermiticity, i.e. the case when hm^y— ►oo(l ~ ''") > Oj t < 1 independent of 
A^. In this regime the averaged density of eigenvalues N~^Ri{Z) is asymptotically 
zero outside the elhpse [ReZ/(l + r)]^ + [ImZ/(l — r)]^ < 1. Inside the ellipse 
lim7v-+oo -^^^^i(^) — [''■(1 — r^)]"^. This sets a microscopic scale on which the 
averaged number of eigenvalues in any domain of unit area remains finite when iV — > 
00. Remarkably, the r-dependence is essentially trivial on this scale: the statistical 
properties of eigenvalues are described by Rn{zi, . . . , z„) = 7V~"i?„(--^, . . . , -^) and 



lim i?„(zi, . . . , z„) = 

N~*oo 



1 



■^?=i'^^''det 



\j,k=l,...,n ■ (16) 



_7r(l-r2)_ 

This limiting relation can be inferred 150! from Mehler's formula for the Hermite 
polynomials 68 . After the trivial additional rescaling z z VI -7-2 the expression 
on the right-hand side in Eq. (|16|) becomes identical to that found by Ginibre |69) . 

Now we move on to the regime of weak non-Hermiticity. We know that in 
this regime new non-trivial correlations occur on the scale: ImZi.2 = 0{1/N), 
ReZi — ReZ2 = 0{1/N). Correspondingly, we introduce new variables x,yi,y2,i^ 
in such a way that: x = Re(Zi -I- 2'2)/2, yi^2 = Nlm{Zi,2), ^ = iVRe(Zi — Z2), and 
consider them finite when performing the limit N —> 00. 

Substituting Eq. I|12|) into Eq. I|15|) and using the above definitions we can explicitly 
perform the limit N 00, taking into account that limjv_,oo -^(1 — ''-)= The 
details of the procedure are given elsewhere 50. . In this regime 

,• K- f , + uj/2-iy2\ / + ix{yi - y2) \ 



1 /•^''-(^) du 



X — / exp 

™ J--nu,,(x) V27r 



? 9 

a u 



u{yi + y2) + iww 



with Vsc{X) = ^\/4 — X"^ standing for the Wigner semicircular density of real 
eigenvalues of the Hermitian part H of the matrices J. 

The kernel Kn given by Eq. p7|l determines all the properties of complex 
eigenvalues in the regime of weak non-Hermiticity. For instance, the mean value 
of the density p{Z) = X^iLi 5^'^\Z — Zi) of complex eigenvalues Z — X -\- iY is given 
by {p{Z)) = Kn{Z, Z*), i.e. by putting yi = y2 and = in Eq.|(T7I) (cf. 03 found 
by the supcrsymmetry approach described in subsection 2.3). 

One of the most informative statistical measures of the spectral correlations is 
the 'connected' part of the two-point correlation function of eigenvalue densities: 

(p(Zi)p(Z2)),- (p(Zi))(5(2)(^^_^2)-3^2(^l,^2), (18) 

In particular, it determines the variance S2(_D) (n(_D)2) — {n{D))'^ of the number 
n = Jjj d^Zp{Z) of complex eigenvalues in any domain D in the complex plane as: 

j:2{d)^ [ d^z, [ d^Z2[{p{z,)p{Z2)) - {p{z,)){p{Z2))] (19) 

J D J D 

= f d^z {p{z)) ~ [ d^Zi f d^Z2 y2{Zi,Z2) 

JD JD JD 

Comparing Ea. (|18|l with the definitions Eqs. H13I15|I we see that the cluster 
function 3^2(-^i,^2) is expressed in terms of the kernel Kj^ as 3^2(-Z^i,^2) = 

\KN{ZuZ*2)t 
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It is evident that in the hmit of weak non-Hermiticity the kernel K^si depends on 
X only via the semicircular density v^ciX)- Thus, it does not change with X on the 
local scale comparable with the mean spacing along the real axis A ~ 1/A^. 

The cluster function is given by the following explicit expression: 



'^''(^^ du 



"(yi + V2) + iuuj 



2 

(20) 



The parameter a = Tiv{X)a controls the deviation from Hermiticity §. When 

a ^ the cluster function tends to GUE form 3^2 (i^, 2/1, 2/2) — ^^(2/1)^(2/2)^^^^—^^^^^^- 
In the opposite case a 3> 1 we expect it to match with the regime of strong non- 
Hermiticity emerging for a ~ N^/"^. To verify this we notice that for large the 
integral in Eq. H20() can be evaluated by the saddle-point method. The saddle-point in 
that case is = — (2/1+2/2 ~ i^) I and to the leading order the result is non- vanishing 
as long as |Reus| < ■kVsc{X') that is \y\ -I- 2/2I < t^Vsc{X^c? . Assuming that this is 
the case one performs the Gaussian integration around Ug and finds (in the original 
variables ^1,^2) the expression equivalent (up to a trivial rescaling) to that found 
by Ginibre 3^2(^1 ,^2) = (A^V7ra2)2 exp{-7V2|Zi - Zap/a^}. Remembering 

VsciX) — 4 — and o? — Nv^ and performing the same procedure for 
the mean density p{Z) — K{Z,Z*) one finds that the condition for the saddle- 
point discussed above just ensures that both Zi and Z2 are inside an elliptic blob 
(Re Z)"^ /A + (2 Im Z)'^ /v'^ = 1 in the complex plane filled in uniformly with eigenvalues 
with the density p{Z) = N/{ttv'^), in agreement with [7D]. This is basically the case 
of "strong" non-Hermiticity, see the discussion preceeding Ea. l(TB|l . and the "elliptic" 
law is believed to be the common feature for strongly non-Hermitian matrices with 
independent entries [HKIITT] . 

One can further calculate the Fourier transform of the cluster function over its 
arguments w, 2/1, 2/2 and find an explicit expression for the spectral form-factor. It 
allows one to determine the variance £2 of a number of eigenvalues in any domain D 
of the complex plane. In this way one traces, in particular, a tendency for gradual 
decorrelation of real parts of the eigenvalues with growing non-Hermiticity. Similarly, 
one can investigate various regimes of eigenvalue repulsion in the complex plane by 
looking at small-distance behaviour of the nearest neighbour distance distribution. All 
these calculations are presented and thoroughly discussed in [511] . 

Further calculations in the framework of the present model were done recently by 
Akemann [321 • Motivated by QCD applications the author first suggested to consider 
the matrix integral 

r r 1 

Z ({m/}) = / dJdJ^ V{J, J^) W Idet (J - imjl) det (J^ -I- imjl) (21) 

J f=i 

with V{J,J^) being the joint probability density function of the random matrices 
J, see Eq.©. This integral is a model for the Euclidean QCD partition function 
in 3 dimensions at nonzero chemical potential, when the chiral structure of the 
Dirac operator is irrelevant. The latter fact gives some justification of the purely 
phenomenological model based on Eq.|5lJ- On the other hand, the matrix integral 
is interesting by itself and its calculation is very intimately related to calculating the 
eigenvalue correlation function Ea. lll|l . 

§ In our earlier Letter |50| we used the definition of the parameter a different by a factor of 2 from 
the present one. 
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Indeed, having in mind the decomposition {Zi, Zpj) ~ ({z}, {^}) as in Eg. ljllll 
and noticing that the Vandermonde determinant A({Z}) can be rewritten as : 



n N — n 



A({z})=A({z})A({a)n n(^''-^:'-^ 

i=i j=i 



it is easy to see that the correlation function Rn{{z}) as defined in Ea. (|ll|l can be put 
in the form: 

2 

TL TL iV — TL 

.1=1 J 1=1 j = l 



CN-n(T) 



(22) 



where VN-nii^}) is just the density of joint probabiUty of eigenvalues ^i, S,N-n of a 
random Gaussian non-Hermitian matrix Jn-h = Hi+ivA of the size (N—n) x (N—n). 
But the integral in the last equation is just the expectation value of the product of 
determinants: 




Za-JN-n) det (z*1-JI 



'N- 



Jn- 



Thus, we see that the problem of evaluating the multipoint correlation function 
Rn{{z}) of complex eigenvalues is, in fact, equivalent to calculating the RMT partition 
function for QCD3, Ea. (|21|l . Akemann|42| used this observation to express the 
partition function Z in terms of the i?„({z}), and in this way exploited the orthogonal 
polynomial representation, sec Eas. 1141151) . 

Here we would like to note that by reversing this kind of reasoning one could 
get access to the multipoint correlations of complex eigenvalues if one were able to 
provide an independent technique of evaluating Z. This was realized by Fyodorov and 
Khoruzhenko |SJ who implemented such an idea for a more general class of ensembles. 
We will see below, that a reduction similar to that discussed above holds beyond the 
Gaussian case, provided one deals with the regime of weak non-Hermiticity. Moreover, 
it will be shown how to perform the ensemble average of the product of characteristic 
polynomials for an arbitrary /ised non-Hermitian deformation by mapping the analogue 
of partition function Ea. l|21|) to a fermionic version of a nonlinear cr— model. 

Let us discuss briefly other types of (weakly) non-Hcrmitian deformations based 
on the Gaussian case and studied by various authors. Recently Akemann |43| 
considered a model with the joint probability density of variables Zi defined by 

Vt\{Z}) cc |A({zf )p |Zi . . . cxp J J2 [l^^l' - i^^i + )] }(23) 

The new factors , a > — 1 as well as squared variables inside the Vandermonde 

determinant were introduced in an ad hoc way to make this ensemble a candidate 
for chiral extension of the Gaussian case, provided one interpretes the variables Z 
as complex Dirac eigenvalues. In such an interpretation the parameter a can be 
related to the number of massless quark flavours Np, whereas the chemical potential 
is phenomenologically identified with ^Vl — t'^- The advantage of the model is 
that it still can be solved by the orthogonal polynomial method, with (properly 

normalized) Laguerre polynomials of complex argument L^,"-* ^-^^f-^ replacing the 
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Hermite polynomials. In particular, elegant explicit expressions for the limiting form 
of the kernel K^°'\Z^,Z*) were derived both in the strong non-Hermiticity and the 
weak non-Hermiticity regimes 03]. On the other hand, the main drawback of the 
model is its ad-hoc construction in terms of the variables Zi. In particular, it is 
completely unclear how to construct the corresponding random matrix model for the 
Dirac operator which could have the starting formula (|23|l as the joint probability 
density of its eigenvalues. The hope is, nevertheless, that the universality will again 
make the eigenvalue correlations insensitive to the particular choice of the underlying 
model H". However, we believe that the full analysis of the random matrix model 
Eq.© in the regime of weak non-Hermiticity still remains an important unsolved 
problem. 

In the spirit of the methods discussed in the present section one can efficiently 
study Gaussian ensembles corresponding to the weakly non-Hermitian quaternion 
random matrices. Matrices with such (symplectic) symmetry are also important for 
QCD applications and that fact inspired quite a few attempts of understanding 
their properties analytically |72[l73l[7^l75j . The most complete solution was presented 
by Kanzieper j75| who used the skew orthogonal polynomials. 

Recently another interesting extension of the model Eq.® was suggested by 
Garcia-Garcia et al.jZHI- The authors introduced an ensemble that interpolates 
between the GUE, the Ginibre ensemble of strongly non-Hermitian matrices and 
the Hermitian Poisson ensemble with uncorrelated real eigenvalues. The physical 
motivation behind such a generalization comes from the wish to understand better 
an interplay between the effects of non-Hermiticity and the effects of Anderson 
localization which drive the statistics of real eigenvalues towards an uncorrelated 
Poisson spectrum. It is claimed that the model can adequately describe critical 
spectral statistics of open disordered systems close to the Anderson transition. 
The authors managed to provide a closed analytical expression for the eigenvalue 
correlations in various regimes and concluded that critical statistics is not modified by 
weak non-Hermiticity effects. They further put forward an expectation that opening 
a critical system does not affect the multifractal dimension of the wavefunction. 
Whether it is indeed the case remains to be verified in numerical simulations. 

2.2. Finite-rank deformations of Hermitian matrices: the case of chaotic scattering 

As it was discussed in the introduction, the poles of the S'-matrix (resonances) are just 
the complex eigenvalues of an effective random matrix Hamiltonian Tie// = H — iT. 
The requirement of the S'-matrix unitarity and causality restricts T to be of the rank 
M, characterized by M positive eigenvalues 7c>0,c=l,..., Af , the rest N — M 
eigenvalues being identically zero. To satisfy the condition of weak non-Hermiticity 
for the case of a few open channels: N M ^ 1 means just to consider these 
eigenvalues to be TV- independent: 7c < oo. The condition ensures the width Tk of a 
typical resonance to be comparable with the mean separation A between neighbouring 
resonances along the real axis. 

Despite quite substantial efforts [HI IT^ ITTl IK^ IHT] non-perturbative results on 
resonance statistics for few-channel scattering were very restricted for a long time. The 
knowledge amounted mainly to (i) the joint probability density of all resonances for 
the system with a single open channel and Gaussian-distributed transition amplitudes 
by Sokolov et al. W (12! ^^"^ (ii) ^^'^ mean density of S'-matrix poles for arbitrary 
M <^ N [111 I52| . A systematic analytical approach to the statistical properties of 
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resonances for systems with broken time-reversal invariance was suggested in j51| . In 
what follows we discuss that approach in more detail. 

For the sake of notational convenience we have changed the sign in front of F ||. 
This trivially amounts to changing the signs of imaginary parts of all eigenvalues 
in the resulting expressions. From now on we consider an ensemble of random 
N X N complex matrices J = H + iV, where H is N x N matrix taken from a 
Gaussian Unitary Ensemble (GUE) of Hermitian matrices with the probability density 
V{H) cx exp (—-J trH^), H = W . As for the matrix F = F'^, we consider it to be a 
fixed nonnegative one: F > 0. 

By analogy with the complex numbers, 

H =l(J + J^) =RcJ; F= -^(J- J^) =IniJ 
2 2i 

with the "only" difference that Re J and ImJ do not commute. Then the probability 
density function in our ensemble of random matrices J can be written in the form 

7'(J)oce-*Tr(Rej)^^(p„j^^j)^ (24) 

with a matrix (5-function for Hermitian matrices. We do not specify the multiplicative 
constant in Ea. (|24(l . It can be found from the normalization condition. Similarly, and 
by the same reason, we will systematically disregard multiplicative constants when 
dealing with probability densities and correlation functions. 

Eq. 1)24(1 can be used to obtain the density of joint distribution of eigenvalues 
by integrating V{J) over the degrees of freedom that are complementary to the 
eigenvalues of J. This again can be done following Dyson's method (see preceding 
section) . To perform the integration over upper triangular matrices R it is technically 
convenient to use a Fourier-integral representation for the (5— function in Eq. H24|l . This 
reduces the corresponding integral to a Gaussian one and after quite straightforward 
algebraic manipulations the resulting expression is 

Pj,{Z) oc e-f I^^Trz=-f Trr^|A(Z)pQM(ImZ), (25) 
where Qm{^''^Z) is the remaining integral 

N 

Qm (IniZ) = / [dU] J] 5 {Imzi - {U^rU)a) , (26) 
1=1 

over the unitary group U{N), [dU] being the Haar measure. It is clear that we can 
consider F to be diagonal, and the integral Ea. (|26|l . has an obvious interpretation of 
the joint probability density of diagonal entries for a matrix with prescribed set of 
eigenvalues F, randomly rotated by a unitary matrix U. Recently this interpretion 
was helpful in applications to the problem of Wigner time-delay distribution I 77| . 

To proceed further we need to integrate over U. Again it is convenient to use the 
Fourier-integral representation for the (5— functions in Eq. H26|) : 

QMilmZ) = 1^ ,^lr.TrKzj [^^^j ^-.TrKU^TU ^ (27) 

where the first integration is over all real diagonal matrices K of dimension N , dK 
being dki . . . dk^. 

II in the rest of this subsection we will denote matrices without "hat" 
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When the eigenvalues 71, ...,7Ar of F are all distinct the integration over U{N) 
can be performed using the famous Itzykson-Zuber-Harish-Chandra (IZHC) formula 
EHl: 



det [. 



where {K} = diag(-ft'i, . . . , Kn) and {7} — diag(7i, . . . , jn) and we have chosen the 
total volume of U{N) equal to 1. We, however, are mostly interested in the case when 
r has a small rank Af ^ N, i.e. it has only M nonzero eigenvalues which we denote by 
71, ... , 7m, the rest N — M being zero. This limit of highly degenerate eigenvalues is 
difficult to perform in the original IZHC formula. A way to circumvent this difficulty 
is discussed in the appendix Al where it is shown that: 

/J i- M—N r N M ^ 

[dU] e-Trxt/trt/ ^ m ^^'^ / dAA(A)e-*5:Jl,7cA. TT TT 1_ (39) 
^(7) J - kj 

with A = diag(Ai, . . . , Am)- Here it is meant that the integration contour is chosen in 
such a way that Im Xc — 0^ . 

To perform the integration over X is a simple task with the use of the identities: 

Ac - fcj ~ ^ Ac - ^^^^ Ac - A,, ' 2^ Ac - A: ~ 1 otherwise 

Taking into account that due to the specific structure of the matrices J their 
eigenvalues lie in the upper part of the complex plane we assume in all formulae 
below that Imzj > for all j. This finally results in: 

QM(ImZ)cxM! ^"7 / dA e^^^- ^=^=A(A) J] E fT^Tx TT ' (^O) 

^(7) f^^ ^ Y\.s{^q) y\ - A.) 

The pair of equations (|25I30|I provides an explicit representation for the joint 
probability density of N resonances Zj in the complex plane. For the particular 
case of a non-Hermitian rank-one perturbation of GUE matrices (single-channel case 
M = 1) the density of joint probability of all complex eigenvalues was found earlier by 
Stockmann and Seba for a closely related model of random couplings W j^. Their 
result follows from our Ea. (|25l30|l by noticing (i) that for M — 1 the A— integration is 
trivial to perform, yielding 

-Af7^/2 N 

7'^(^)|M=i(xe-fReTrz^|A(Z)p^^^J(7-^Imz,) (31) 

^ i=i 

and (ii) in their case the only nonzero eigenvalue 7 > is a random variable with the 
probabihty density 7^(7) esc 7^"^ exp — [A^7/7o]. 

Our main goal is to use equations H25I30|I as a basis for calculating the n-eigenvalue 
correlation functions i?„({z}) as defined in Ea. HllI) . In what follows we will calculate 
Rn{{z}) for arbitrary fixed n, M in the limit N — > 00, but it may be instructive for the 
reader first to verify all the steps on the simplest case Af = 1, starting from the Eq.JSTJ 
and following [ST]- On the first stage we will replace the integration over f in Hll|l by 
averaging over the ensemble of non-Hermitian random matrices Jn-u = -ffw-n + *r, 
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with Hj^^n being a GUE matrix of the reduced size {N — n) x {N — n) (cf. Eq.||22)). 
In this way one finds the foUowing representation for the correlation functions: 

where we introduced the notation gc = 5 (7c + Ic^) ''^^'^ denoted: 

C^iM) - (n I det [z, - JN-nil)] n (33) 

/ GUE 

The corresponding derivation is presented in appendix B and is vaHd in the Hmit 
n,M < ^ 00. 

Thus, the problem amounts to evaluation of the correlation function of the 
determinants in Eq.(j2SJ- To proceed we first write each of the determinants as a 
Gaussian integral over a set of Grassmann variables XiX^'- 

n I det [zp ~ Jjv-n(7)] r = / dxdx^e-(x'^^"^+^'H"-">^+*'^'[^®^^"]^) . (34) 
p=i 

When this is done, the GUE average over Hm-u becomes trivial and yields the 
terms quartic with respect to the Grassmannian vectors. These terms can be further 
traded for an auxilliary integration over a Hermitian matrix S of the size 2n x 2n 
(the so-called Hubbard-Stratonovich transformation). Then the integration over the 
Grassmann fields becomes again Gaussian and is trivially performed yielding again 
a determinant. Introducing diagonal matrices defined as: Z2„ = diag({z}, {z}''') and 
L2n = diag(l„,— 1„) and using the finite-rank property of F one easily brings the 
result to the following form: 

» M 
C^{{z}) = / [d^]e-(^-")Tr[i5^-in(Z.„-..s)] ^^^^ ^ Z7,L2„(Z2„ - iS)-^] .(35) 

c=l 

We are interested in evaluating the above integrals in the limit — > 00. Let us 
now recall that nontrivial eigenvalue correlations are expected to occur on such a scale 
when the eigenvalues are separated by distances comparable with the mean eigenvalue 
separation for GUE matrices H (see preceding section), the latter distance being of 
the order {N — n)~^ for our choice of V{H). Accordingly, it is convenient to separate 
the "center of mass" coordinate x — - Y^" , Rez, so that Zj ~ x + .f' , where both 
the real and imaginary parts of Zj are of the order of unity in the limit when N 00 
and M is fixed. In this limit i?„({z}) is effectively a function of {z} {x is fixed) which 
we are going to calculate. To this end, we expand the logarithm in the exponent of 
Eq.jnS) as: 

{N - n)Trln(Z2„ - iS) « (TV - n)Trln(2;l2„ - iS) + Tr [Z2„(a:l2„ - tS)-^] + ...(36) 

keeping only terms up to the order 0(1) when N —^ 00. With the same precision we 
set Z2„ = xl2n in the pre-exponential determinant factor of the integrand, Ea. (|35ll . 

The next step towards evaluating the integral in H35|l is to diagonalize the 
Hermitian matrix 5" as: S — C/2nSJ7^^, where U2n G U{2n) and E — diag(o'i, o'2n)- 
Then, keeping only the terms relevant in the limit N 00, we obtain: 

r ^ -(Af-n)l: [^-ln(x-4(Tfc)] 

C^({z}) = ydSA2(S)e -1 {C{S))ui2n) (37) 
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where 



{C{S)) 



(7(2ri 



= /M[/2„]e-Tr[Z.„(. 



M 



"'1 Y[ det [l2« + z7cL2„(xl2„ - i^)"'] (38) 



c=l 



The form of the integrand in (|37|l suggests exploiting the saddle-point method 
in the integral over a^, k — l,...,2n. Altogether there are 2^" saddle-points: 



= —^{x + iek^J A — x^) each corresponding to a particular choice of 2n signs 
±1. It is easy to understand, however, that not all these saddle points are 



equally important in the limit of large N . The leading order contribution comes from 
integration around those saddle-points where exactly n parameters equal 1 (the rest 
being equal -1). All other choices can be neglected as they lead to lower order terms. 
This is because of the presence of the Vandermonde determinant in the integrand: 



ki<k-2 



ki<k2 ^ ^ 



where we denoted by Sak fluctuations around a given saddle-point to be treated in 
the Gaussian approximation. It is now evident, that we should require that the 
minimum number of the differences ek^ — Cfea: /ci < fc2 = 1, 2, n vanish. Otherwise, 
the Vandermonde determinant is proportional to higher powers of Sak, each power 
producing a small factor {N — n)~^ when integrating around the saddle-point in 
the Gaussian approximation. This requirement fixes the choice of the saddle-points 
described above. 

At the same time, all relevant saddle-points are equivalent (because they can be 
transformed one to the other by an element of the U{2n) group) and produce the 
same contribution. Therefore, we can take one of them, suppressing the multiplicative 
constants as usual. Our choice of S*^*) is S^*^ = — |xl2n + T^v{x)'L2m where the 
symbol i'{x) stands for the semicircular density of real eigenvalues of the matrices , 
^(^) ~ 5^\/4 — x^. Substituting E*^**^ to the integrand it is easy to verify that in the 
limit — > oo we have: 

C7^({z}) cxe^S^=i^°^'q({i}) 

where {z} = (zi, . . . , z„), Zj ~ {N — n)(zj — .x), j = 1, . . . , n, and 

M 



(39) 



.(a;)TrZ2„Q, 



Yl det 



c=l 



7ri/(a;)7cL 2nQ2„ 



(40) 



In (|in|l Q2„ ~ U2n1^2nU2n and the integration is over the coset space U{2n)/U{n)® 
U{n). 

Thus, the problem reduces to evaluation of an integral over a coset space. This 
type of integrals is known in the literature under the name of zero-dimensional 
nonlinear cr— models and the corresponding calculation is a standard one and is 
outlined in |51| . The result is given by: 



C^Mz}) (X 



det 7" 



n\ det 7" 

mm- 



det 



dXr, det 



M 



det 



n M 



3=1 c=l 



dX Y\_i9c + 7Tiy{x)X]e 
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Combining this with Eqs. (|32|) and (|39|l . denoting u = XTrh'{x) and restoring the 
normaUzation we finaUy see that the correlation functions have the foUowing simple 
structure: 

^ Rn(x+^,...,x + ^) =det[K(~z„rk)]h.k=i.,...,n , (42) 



where the kernel K{zj, zl) is given by 



M 



1 p-KV 

KCz,, z;) = -Fi/2(5i)Fi/2(i2) / duJligc + u] 

with F{z) — X^fcli U"^ — (g !!°g^.) ^(Imz). Introducing a "characteristic function" 

M . s 

/rM = ^ln(l + ^j (44) 



dfc exp(-ifcImf-/r(-iA:/2))(45) 



and noticing that 

F(f) = 0(lmz)^=-^- n- r/ 

we can rewrite the kernel Ea. 143(1 in an equivalent form. Denoting zi = iyi—Lu/2; Z2 ~ 
iy2 + uj/2 we find that the kernel is given by : 

K(Zi, Z*) ^ N^K(Si, zn*) = ^ / du e-(«i+!'2)«+i'^«+/r(«) 

/ poo nOQ 

X ( / dfci exp{-jfci?/i - /r(-«fci/2)} / dfc2 exp{-ifc22/2 - /r(-ifc2/2)} 

Using the derived expression we are going to demonstrate that for a system with 
many open channels M 3> g the statistics of resonances eventually becomes Ginibre- 
like. This fact coordinates with existing numerical simulations [161 117L ESI • 

To this end we consider the case of equivalent channels: gi = g2 = ■■■ = gM = g 
and evaluate all the three integrals in Eq. ((46|l by the saddle-point method in the limit 
Af 3> .9 ~ 1. A straightforward calculation gives 



/oo 1 
dki exp{-zfci2/p - fr{-ikj,/2)} « M^/^-^' e^' ^/2^-{2gyj,)^' e-^^y- 
-oo Vp 



for p = 1,2 whereas the integral over u dominated by the saddle-point Ms = 
—g + M/ys ; ys — yi + y2 — ii^ is estimated in the limit A/ — > cx) as 

M^'+^/\gys)-^'e-^'+ay^V2^- 

Vs 

as long as |Reus| < tti/. In the opposite case |Reus| > ttv the integral is zero. 
Combining these expressions and taking the absolute value of the kernel results in: 

\m..Zl)\=N- ^ X(^y,Ty2)K.-r^^-^l- ' ^^»^(47) 
This expression immediately shows us that the mean density of resonances is 

p{Z) = \KiZ,Z*)\=N^—^^—-—^ ; -— -<y< 



47rj/2 47r(ImZ)2 2{g + nv) - 2{g - nv) 
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and zero otherwise. Such a "cloud of resonances with a gap" is in agreement 
with the earlier result for the mean density obtained by Haake et al.^^l in the 
limitM ^ c» ; — > oo in such a way that M/N = fixed. 

Let us now suppose that both Zi and Z2 are within the " cloud" . Then it is 
possible to rewrite Ea. H47|l as: 



\k{z,,z;)\^{p{z,)p{Z2)) 



M + l 

1/2 f-. _ {yi-y2)^ +uj^^ 



Taking into account that actually yi^2 > ^11^9 + Stti/) 1 we see that always 
2/1 + 2/2 S> ^ 1 whereas the relation between y\ — j/2 and lo can be arbitrary. 
Introducing 



P 



we see that 



\k{Z,,Z;)\ = {p{Z,)p{Z2)f^ [l ~ ^^^^[iy, - y2r + LO^]j (48) 

-^{p(Z,)piZ2))'/'e-^M^)\^^-^^\' 

in the limit M — » 00. It is further evident, that the correlations vanish very rapidly 
when the distance \Zi — Z2I between eigenvalues exceeds M~^/^. Thus, we can safely 
put p{Zi) w p{Z2) ~ p ( ^^+^2 ) when M ^ 00 arriving finally to the Ginibre-like 
expression for the kernel: 

\K{Z,,Z;)\^p{Z)e^p-^7Tp{Z)\Zi-Z2\^ (49) 

In contrast to the Ginibre case, here the density p{Z) inside the cloud is not constant 
but is rather position-dependent. It is natural to conjecture that expression Eq.|03 
constitutes the most universal form of the eigenvalue statistics for the case of strongly 
non-Hermitian random matrices. This conjecture is also supported with existing 
numerics, e.g. for a QCD-inspired non-Hermitian model |40| . 

Having at our disposal the correlation functions, we can easily repeat the 
calculations done in the previous sections and find such statistical characteristics as 
spectral form- factor and the number variance |^ . In particular, for the case of 
equivalent channels gc — g, c ~ 1, M the variance of the number of resonances in a 
strip < ReZ < L = L^A; ~oo < Im Z < 00 is given by: 



1 



(.9 + vf 



(50) 



dv 

'-(i-fc) 

where we have used 'kv{Q) = 1. Let us discuss for simplicity its typical features for the 
simplest case of a single open channel M = 1. We are actually interested in deviations 
of the number variance from its value known for Hermitian GUE matrices. One finds: 

5E2(i.) = S2(L,) - S^^^(ix) = A / ^^sin^ [TTkL^) ^ \7^.., (51) 

As usual, we are interested in the behaviour of the number variance for L^, 3> 1, where 
Ti2^^{Lx) grows logarithmically as: [Lx) ~ InL^. -f const -f ... For any g > 1 

we find that the difference (5S(Lj.) tends to a constant value In [5^/(2/^ — 1)]. This 
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fact means that the curves ^^'^^(Lj.) in the asymptotic regime L^: ^ 1 are just shifted 
upwards by a constant amount. The shift is the larger the more "open" is the system, 
i.e. the closer is the value of g to unity. The so-called " perfect coupling" case g = 1 is 
known to be in many respects specific . Physically, it describes the situation when 
direct scattering is completely absent. On a level of number variance such a specificity 
is reflected in changing the coefficient in front of the leading logarithmic growth from 
the GUE value I/tt^ to a new value 3/(27r^). 

Qualitatively, the same picture holds for M > 1: for a fixed value of the strip 
widths Lx the number variance is the larger, the larger is M and the closer are the 
coupling constants g to unity. This fact is just a fingerprint of a gradual " decorrelation" 
of resonance positions along the real axis. As to the small-distance behaviour of the 
nearest-neighbour distance distribution p{Zo,S ^ A) in the complex plane , the 
leading term for S ^ turns out to be always cubic: p{Zo,S <C A) oc 5*^, as long 
as the system is open, in agreement with the existing numerical data ^| I18L . 
However, for few open channels and 5 3> 1 the interval of cubic behaviour turns 
out to be parametrically small: S ^ g~^^- For the observation points Zq taken at 
regions with asymptotically low density of resonances a crossover to an anomalous 
S*^/^ behaviour might be observable, cf. pT)] . 

In the end of the present section we discuss briefly an intimate relationship 
which actually exists between the kernel Eq. (|43|l derived for a flnite-rank perturbation 
and its Gaussian counterpart, Ea. lfTTIl . For this we notice, that the kernel Eq.lQJJ 
derived under an assumption of positivity of each 7c ( hence gc), but in fact its 
form Ea.H46|l describes the statistics of complex eigenvalues for arbitrary finite-rank 
perturbations, irrespective of the sign of a particular eigenvalue 7c. This fact can 
be proven by a straightforward modification of the technique outlined in the present 
section. Moreover, a little thought makes it natural to expect that its validity is not 
restricted even by finite-rank perturbations, but rather comprises the whole class of 
almost- Hermitian matrices, characterized by A^TrF^ cx triJ^ in the limit N ^ 00. 
As a confirmation of this conjecture, let us look from this point of view on a typical 
Gaussian deformation. For that case a typical eigenvalue of F is 7c ~ N~^/'^ 1 and 
therefore ~ 27c, c = 1, ...,iV. This fact allows us to expand the "characteristic 
function" in Eg. lB^ as: 

/r(u)k^oo ~ 2utTT - 2u^TtT^ + .... (52) 
It is an easy task to show that for a Gaussian F we have typically TrF ^ N~^/^ ^ 
TrF^ ~ 0(1). Therefore one can put safely: fr{u) ~ — 2w^TrF^ which renders the 
/c— integrals in Ea. l|4H|l to be Gaussian ones. Identifying: TrF^ = and performing 
the Gaussian integrals we see that the resulting kernel coincides with that given in 
Eg. !]!?!) , up to an overall phase factor e'^/^f^^^"^!) which in any way does not play any 
role in calculating the correlation functions of eigenvalues. 

2. 3. Density of complex eigenvalues via supersymmetry approach 

Unfortunately, it is not known at the moment how to extend the methods described 
above to the case of non-Hermitian deformations of real symmetric matrices. We 
consider this issue as one of the most challenging problems for future investigations. 
At the same time the mean density of complex eigenvalues turns out to be accessible 
via a technique known as Efetov supersymmetry approach and is shortly discussed 
below. In fact, the method provides access to non-Hermitian deformations of matrices 
interpolating between real symmetric and complex Hermitian, as discussed in |50|. 
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As before we decompose any N x N matrix J into a sum of its Hermitian and 
skew-Hermitian parts and consider an ensemble of random N x N complex matrices 
J = Hi + ivH2 where Hp; p = 1,2 are both Hermitian: Hj^ = Hp. The parameter v 
is used to control the degree of non-Hermiticity. In turn, complex Hermitian matrices 
Hp can always be represented as Hi — Si +iuAi and H2 — S2 + iwA2, where Sp = S'J 
is a real symmetric matrix, and Ap = —Ap is a real antisymmetric one. From this 
point of view the parameters u, w control the degree of being non-symmetric. 

In the first part of this subsection we consider the matrices 51,52,^1,712 to be 
mutually statistically independent, with i.i.d. entries normalized in such a way that 

hm ^TrSl = hm ^TrApA^ = 1 (53) 

Such a normalization ensures that for any fixed value of the parameter u > statistics 
of real eigenvalues of the Hermitian matrix H = S + iuA in the limit — > 00 is 
identical (up to a trivial rescaling) to that of u = 1, the latter case being standard 
GUE. On the other hand, for w = real eigenvalues of the real symmetric matrix 
S follow a different pattern of the Gaussian Orthogonal Ensemble (GOE). The non- 
trivial crossover between GUE and GOE types of statistical behaviour happens on a 
scale u oc [T]. 

Similar arguments show 07| |5nj , that the most interesting behaviour of complex 
eigenvalues of non-Hermitian matrices should be expected for the parameter v being 
scaled in a similar way: v cx iV~^/^. It is just the regime of weak non-Hermiticity 
which we are interested in. Under these conditions a non-Hermitian matrix J still 
"remembers" the statistics of its Hermitian part Hi. As will be clear afterwards, the 
parameter w should be kept of the order of unity in order to influence the statistics 
of the complex eigenvalues. 

Correspondingly, we scale the parameters as \. 

(54) 



2\/]v' 2Vn 
and consider a, 4>, w fixed of the order 0(1) when A^ 00. 

For calculating the mean density of complex eigenvalues Zk = Xk + iYk, k ~ 
1, 2, ...,N defined as 

N N 

p{Z) = S^^^ {Z -Zk) = J2 - ^k)S{Y - Ffc) = p{X, Y) (55) 

k=l k=l 

one relates it to the " charge potential" : 



$(X, Y,K) = j- \nDet[{Z - n,jf){Z - HeffY + n^] 

in view of the relation: p{X,Y) = lim^^o 5^ $(A, F, k), where stands for the 
two-dimensional Laplacian, see appendix A of Technically, it is convenient to 

introduce the generating function (cf.JTj) 

2 ^ De^ [{Z-n,ff){Z-neffV + n^] 

Bet[{Zt~neff)iZt-neff)^ ^ ' 

in terms of which 

\ d d 

p(Z) = lim - — lim —2 (57) 

' TT K^o dZ* z,-^z dZb 

^ In the Letter I5UI there is a misprint in the definition of the parameter a. 
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Assuming the matrix elements of all involved matrices to be statistically independent, 
up to symmetry constraints, one performs the ensemble averaging by the standard 
trick of representing the ratio of the two determinants in Eq. H5t)|l in terms of Gaussian 
superintegrals. After a set of manipulations exposed in ,50 | one arrives at the following 
expression: 

y)) = J d^l{Q) Str (af )q) Str (a.g) exp -S{Q) (58) 

2 2 ?i2 2 2 2 

5(Q) = -^yStr(<T,Q)-^Str(a,g) + -Str (f2Q) - ^Str (^q) (59) 

where we introduced the scaled imaginary parts y = ttv{X)NY , with I'iX) being 
the mean (semicircular) eigenvalue density. We also used the notations: = 
{TTL'{X)a) , b'^ — {nu{X)(j)) , = {TTi'{X)aw) . Here the integration goes over 
the set of 8 X 8 supermatrices Q satisfying the constraint — —1. The symbols 
Str, Sdet stand here for the graded trace and the graded determinant, correspondingly. 

Properties of these matrices and the integration measure dfj,{Q) can be found in 
|29| . The matrices ai^'' , ct^ and f2 in Ea. (|58|l are some fixed block-diagonal 8x8 
supermatrices, and can be found in ^EU^. The whole expression is just a general a— 
model representation of the mean density of complex eigenvalues in the regime of weak 
non-Hermiticity. We expect it to be universal and applicable even beyond the case 
of matrices with independent entries, after replacement of the semicircular density 
v{X) with actual density of real eigenvalues of the Hermitian part. The expression 
is parametrised by a, 5 and c. The parameters controll the degree of non-Hermiticity 
(a) , and symmetry properties of the Hermitian (b) and non-Hermitian (c) parts of the 
ensemble. 

Still, in order to obtain an explicit expression for the density of complex 
eigenvalues one has to evaluate the integral over the set of supermatrices Q. In general, 
it is an elaborate task due to the complexity of that manifold. At the present moment 
such an evaluation was successfully performed for three "pure" cases: those of almost- 
Hermitian matrices |47| . real almost-symmetric matrices |45| and complex symmetric 
matrices The first case (which is technically the simplest one) corresponds to 

(/) — > cxD, that is 6 — !■ oo. Under this condition only that part of the matrix Q which 
commutes with T2 provides a nonvanishing contribution and the resulting integral j47| 
proves to be equivalent to the "diagonal" part K{Z,Z*) of the kernel Ea. (|17|l . as 
expected. 

The second nontrivial case for which the result is known explicitly is due to 
Efetov 02]. This is the case of real matrices with "weak asymmetry". In the present 
notations this case corresponds to the limit — > 0; w — > cx) in such a way that the 
product (j)w = c is kept fixed. The density of complex eigenvalues turns out to be 
given by: 

px{y)^6{y) f exp (-52^2/2) (60) 
Jo 

+2y^-y^ dwexp^-^^^^^ dttsmh{2t\y\)exp{-d!^t^/2), 

The first term in this expression shows that everywhere in the regime of "weak 
asymmetry" c < oo a finite fraction of eigenvalues remains on the real axis. Such 
a behaviour is qualitatively different from that typical for the case of "weak non- 
Hermiticity" a < oo, where eigenvalues acquire a nonzero imaginary part with 
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probability one. In the limit c >> 1 the portion of real eigenvalues behaves like 
c~^. Remembering the normalization of the parameter v, Ea. H53|l . it is easy to see 
that for the case of v = 0{1) the number of real eigenvalues should scale as y/N. The 
fact that of the order of N'^/'^ eigenvalues of strongly asymmetric real matrices stay 
real was first found numerically by Sommers et al. |70[ I79| , and proven by Edelman 

At last, the case of complex symmetric matrices (6 = c = in the present 
notation) turns out to be most involved technically. The result of evaluating the 
integral Ea. (|55|) over the coset space turns out to be: 

where 

^^--'^^"^'^(T^l^yTTi- (62) 

Actually, the integrals over Ai , A2 as they stand in Eq. should be understood 
in the sense of a principal value. To be precise, the correct expression is equal to 
the half-sum of two integrals with integration contours encircling the singular point 
Al = — iA; A2 = —i\ from above and from below. 

Of course, the supersymmetry approach is not restricted to anti-Hermitian 
deformations iT whose entries are statistically independent as assumed above, but 
can be used as well for fixed finite-rank deformations. This allows one to calculate the 
density of resonance poles in the complex energy plane. In fact, such a calculation 
was the first application of the supersymmetry approach to non-Hermitian problems 
[T5) and paved the way to all later developments in j37| 23 EH ■ 

One finds the following cr— model representation for the distribution of scaled 
resonance widths y = ^^'^^ < (measured in units of the local mean level spacing A 
of the closed system) for the resonances whose positions are within a narrow window 
around the central point X = of the spectrum. 

{Px{y)) = ^ y d^l{Q) Str (af )q) Str (a,g) exp ^yStr (a.g) (63) 

M r . n 

X II Sdet-1/* l_^|g,a,| 

a=l L ^'^ 

Here |(5, (7r| = Q^t + (^tQ stands for the anticommutator. 

Performing the integration over the manifold of the supermatrices Q (which is of 
different form for different symmetry classes) one finds for the case of broken invariance 
| 15| the distribution coinciding with that following from the kernel Ea. (|43l) . whereas 
for the case of preserved invariance the distribution is given by Ea. (|61|l . but with the 
function F(A) replaced with the following expression (S2j: 

1 

^'^'^ = mii(.c-A)V.(l-A^)V. (64) 

Such a replacement looks completely natural from the point of view of the 
correspondence pointed out in the end of preceding section, see Ea. (l52ll . 

Finally, we would like to mention a finite-rank generalization of Efetov's result, 
Eg. HftUI) as performed in j61| . Namely, we consider the ensemble of real asymmetric 
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matrices of the form: H + A, with H being a real symmetric (taken from GOE) and 
A being a fixed real antisymmetric: Aij = —Aji of a finite rank M . As is well-known, 
a general antisymmetric matrix A of even dimension can be brought by an orthogonal 
rotation to the following block-diagonal form: diag(^i, ....,^jv), with each block Ai 

being 2x2 matrix of the form Ai ~ ( ^ 'n ) • Invoking the argument of the 

rotational invariance of GOE matrices, it is enough to consider deformations to be of 
such block-diagonal form. We restrict our attention to the finite-rank case, when all 
jii = for i > M. Introducing the scaled variable y = Tri'{X)NlmZ we find for the 
rescaled density px{y) = {p{Z))W''^'^{^)]^^ the following expression 

px{y)--5{y) I -h — / dttsinh(|y|t)e5[^^('^^*)+^^(-'^'^*)l (65) 





oo /"OO 



27r 



A 



ds I dkexp{~iks [/^(iTrz^fc) + /^(— iTrt^fc)]} 

\y\ "'-oo 2 

where the function fA{z) is given by the same expression as Eg. 1441) . but with 7c 
replaced by pi, i.e. gc replaced by (/i^ -I- l/pi)/2. 

Our interest in the problem was stimulated by the papers claiming 
that expression Ea. H60(l . derived for Gaussian A, describes well the complex 
eigenvalue density even for non-random antisymmetric deformations by matrices 

^ ~ ), with a constant p being of the order of p ^ N^^^'^. This 

fact immediately follows from Ea. H()5|l . Indeed, for M oc N and the typical pi of the 
order of N~^/'^ the function fA{v) can be expanded up to a first non- vanishing order 
such that fAiv) + fA{—v) oc w^TrA^, cf. Ea. l(5^ . The corresponding expression then 
indeed coincides with Ea. (|()n|l . 

3. Eigenvalue Statistics of Subunitary Matrices 
3.1. Truncations of random unitary matrices 

Let us first consider truncated matrices drawn from the circular unitary ensemble 
(CUE). Let U ^ n ) be an iV X matrix from CUE and A a subunitary 



{N — M) X (TV — M) matrix. To make contact to the general case discussed in 
the introduction we note that putting B — and D — corresponds to a matrix 
A = u\/l - f with T, = for i = 1, ...,N - M andT^ = 1 for i = - M + 1, ...,N. 
This matrix has the same nonzero eigenvalues as the subblock A of the matrix U 
and we determine below the joint density of those eigenvalues inside the complex unit 
circle. 

The joint density of elements of U can be written as 

P{U) cx S {A^A + C^C ~l)S {A^B + D) SiB^^B + D^D - 1) (66) 

with appropriate matrix ^—functions. Integrating out B and D we obtain as joint 
density of elements of A 

P(A) cx /dC (5(AU-|-C^C~ 1) (67) 



with a 2M{N — M) dimensional integration over the complex parameters C. 
Performing these remaining integrations is possible after representing the i5-function 
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in Ea. (|67|l as a Fourier transform over the set of Hermitian matrices J- of the size 
(TV — M) X [N — M). Then integrals over G are essentiahy Gaussian and can be 
immediately performed yielding: 

P{A)oz J dJ' detJ"-^' exp{iTTj^{l- A'fA)} (68) 

The matrix integrals of this type emerge frequently in various context, and most 
recently were calculated in |^|^. For the case M > N/2 the probability P{A) is 
given by a simple, non-singular expression: 

P{A) cx det{l - #A}2*^-^6i{l - AU} (69) 

where the matrix 0— function is unity when its argument is a positive definite matrix 
and is zero otherwise. 

In the opposite, most interesting case M < N/2 the integral still can be calculated 
explicitly |92| but the result is complicated and contains a product of many singular 
((5— functional) contributions. 

We are however mostly interested in finding the joint probability density of 
complex eigenvalues rather than just distribution of the matrix A itself. This can be 
achieved directly from the equation Ea. (|67|l if we again use the Schur decomposition 
A = V{z + R)V^^ extensively exploited by us in Sec. 2. Here the transformation 
V is unitary, z is a diagonal matrix of the complex eigenvalues of A and R is 
strictly upper triangular. A procedure allowing one to integrate out matrices V, 
the {N - M){N - M - l)/2 complex parameters i?y and the M{N - M) complex 
parameters C is described in some detail in J^. The resulting expression is however 
rather simple: 

1...N-M N-M 

p({z})(x n n (70) 

This distribution is very analogous to the Ginibre ensemble and all correlation 
functions can be found readily by the method of orthogonal polynomials. In fact, 
the powers z"^-'^ are already orthogonal inside the imit circle of the complex plane. 
An equivalent method is to consider the joint density P{{z}) as the absolute square 
of a Slater determinant of normalized wave functions 

0„(z)=z"-i«;(|zn/VA^ 

, with Nn standing for a normalization factor and the domain of integration being 
|z| < 1. The kernel, which determines all correlation functions is given by 

N-M 

K{z,,z*)= {z,z*r-^w{\z,\')w{\z,\')/N^. 

n=l 

Further defining x = r"^ = |zp we can, for example, easily write down the averaged 
density of complex eigenvalues at point |z| as 

Pir) . (71) 
^ ' N -M {M \dx) 1-x ^ ' 

There are two important limiting cases for large N: either /j- = M/N fixed or M 
fixed. For fixed /i and — s- oo we find the scaling behaviour: 
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for < 1 — ^ and P{r) — otherwise. The distribution shows a gap near the 
unit circle and quahtatively resembles one obtained for resonances in the chaotic 
scattering problem for large number of channels, see Haake et al.^3|- In this - 
"strongly nonunitary" - limit one can also simplify the cluster function given by 

y(zi,Z2) = |i^(^i,^2*)l': 

z + (5) = (AT - Mfp[zf exp(-7r(iV - M)p(z) (73) 

This is just the Ginibre behaviour Ea. H49l) with the distance i5 rescaled by the local 
mean level distance Xj \J(N — M)p{z) given by (ff^ through p{z) = P(r)/27rr. The 
same can be shown for the nearest neighbour distance distribution obtained by Grobe 
et al. IHOl and applied to a damped chaotic kicked top. 

In the other limit of fixed AI and iV to oo the matrices may be considered as 
weakly nonunitary, and we recover exactly the universal resonance-width distribution 
^1^2 foi' perfect coupling to M channels with y = N{1 — r) 

Similarly the cluster function obtained in this limit can be shown to coincide with the 
one obtained in the previous section for chaotic scattering with a finite number M 
of perfectly coupled channels. The statistics (|74(l has also been found by Kottos and 
Smilansky ,22, for chaotic scattering on graphs and by Gliick et al.,23, for a model 
of crystal electron in the presence of dc and ac fields. In both of these works the 
^-matrix is reduced to the resolvent of a subunitary matrix as is investigated in the 
present subsection. 



3. 2. General finite-rank deviations from unitarity 

In what follows we assume all eigenvalues of the deviation matrix T = 1 — G are 
smaller than unity: < 1, but the resulting expressions turn out to be valid in the 
limiting cases some = 1 as well. 

We again start with the Schur decomposition A = U(Z + R)U^ of the matrix A 
in terms of a unitary matrix U, a diagonal matrix of the eigenvalues Z = {z} and a 
upper triangular matrix R. One can satisfy oneself, that the eigenvalues zi, zat are 
generically not degenerate, provided all < 1. Then, the measure written in terms 
of new variables is given by dA = \ A{{z})\^dRdZdii{U), where the first factor is just 
the squared Vandermonde determinant of eigenvalues Zi and dpiU) is the invariant 
measure on the unitary group. The joint probability density of complex eigenvalues 
is then given by: 

V{{z}) cx |A({z})|2 j dp{U)dR 5 ({Z + R){Z + R)^ - U^Gu'^ . (75) 

The integration over R can again be performed using its triangularity (similar to the 
case of "truncated" matrices considered above). We give here a derivation in more 
detail. Let us consider the columns of the N x N unitary matrices U as A^-component 
(mutually orthogonal) vectors a/, I — 1, N , so that U — (ai, a2, ajv). Then the 
5- functions in Ea. (|75|) imply 

|z,|2+ ^ |i?Hp+alfa, = 1 (76) 

k{<t) 
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and for i < j: 

z:R,, + J2 KrRkj + alTa, = (77) 

k{<t) 

Integration over the complex parameters Rij can easily be done and yields the Jacobian 
l/ni<j l^iP- The remaining integration over the unitary group should be done with 
the product of the diagonal (5-functions as integrand. There the solution i?y of Ea. H77|l 
has to be inserted. We can find that solution employing the hierarchical structure of 
the equations. It can easily be checked, that the solution is given by 

Ri^^ = - 3^4—53 ^ ^Ta., . (78) 



z 



k 



^l- f (ai «) a| + ... + afe_i ® a^^_^ 



One may also show that Eq. H77|l can be extended to z = j such that 



k<i 

- ot 

1 - T(ai ® a{ + ... + a,_i ® al_^) 



= 4 . ^, „ t . . -^^a, . (79) 



This implies 

P({z})(x|A({z})p ^ (80) 

Ilj<j I ^2 I 

X [dix{U)X{5{\z,\''-l + 4—-^ ^ T^ra,). 

J \ 1 - T(ai (8)ai + ... + ai_i (8)aJ_i) 



Finally we use the identity 



t 1 ^ _ det(l-T(ai®a{ + ... + a, ®aT)) 

L — Si - 7 7 J- 3.7 



1 - T'(ai (g) a| + ... + a,_i (g) aj_ J det(l - T(ai » a| + ... + a,_i (g) aj_ J) 

and introducing the projector matrices Pi = diag (1, ...1, 0, 0) after a simple algebra 
we arrive at the following expression: 

P({z})cx|A({z})|2 /d;,(J7)n'5(kir--kd'-det[l-TOP,;7t]) . (81) 
•' 1=1 

The remaining integration over the unitary group can be done, in essence, similar 
to the corresponding procedure for non-Hermitian matrices, see Section 2.2. Actual 
calculations are, however, more involved and we outline the main steps of the full 
solution below extending 

First, introduce x iV matrices of rank I: Qi — a^iga- , so that (jPiU^ = Qi. 

Due to the fact that only M out of N eigenvalues of the matrix T are non-zero, both 
the matrices Qi and the vectors a can be effectively taken to be of the size M (it 
amounts to changing the unspecified normalization constant in Ea. (|81|l 'l. We then 
redefine the matrix T as T = diag(ri, ...,Tm) = t^t, i.e. the matrix of transmission 
coefficients < Tq, < 1. 

Writing down the corresponding constraints in a form of 6— functions, one can 
represent the expression Ea. 17511 in the form: 

. N N 

V{{z}) cx I A({z}) p / n rf'ai S{1 - ^ aj ® a]) 

i=l j=l 
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N 



1=1 



dct 



(82) 



Here ai, i — 1,...,N are AZ-dimensional complex vectors building a, N x M 
subblock of a unitary matrix U and we define (Paj := J^^f^j^ Redai^-') Imdai^^^. 
The unitary constraint is forced by the inserted (5- function. Obviously jz^p — 

det{l-fj:W 



a)/det(l-TE; 



al) < 1 as it should be for complex 



eigenvalues of a contraction. 

We want to integrate over the unitary group which amounts to integrate over all 
complex vectors &{. We increase again the number of integration variables introducing 
as independent variables matrices qi through the dyadic product qj — (g) aj. The 
M X M Hermitian matrix qi contains in general more independent real variables than 
the complex vector ai. Therefore we insert the identities 1 = / dq^ 5{qi — ai (g) a!) 
with integrations over the Hermitian matrices qi and the corresponding J-functions 
for independent elements. Then 182|l takes the form 

„ N N N 

V{{z})^\A{{z})f J l[dq.,Sil-Y,qj)l[ J rf'a<5(9,-a0at) 



N 



\zi\ 



\zi\ 



det 



1=1 



T 



i=l • 
I 



(83) 



Finally, we introduce Qi — 9i ^'^id observe that Eau.f lS^ aquircs the form 



P({z})(x|A({z})p / H'^^Oll'^O^il 




1=1 



zi 



a (Ki a' 



det(l - TQi 



(84) 



where the matrices Qi are considered to be unconstrained N x N Hermitian. We also 
used the orthonormality condition Qn — 1 as well as the convention Qq = 0. 



3.3. Rank-1 deviations from unitarity 

The remaining integration over the vectors ai still requires a quite lengthy calculation. 
The general case will be considered in the next section. A particular case which 
can be solved relatively easy is the case of a rank-one deviation from Hermiticity: 
M — 1, corresponding to a system with only one open channel. This example is quite 
instructive and the exposition below follows ■ Evaluating the integrals Eq. ( I84|l 
quite straightforwardly we arrive at a very simple expression: 

V{{z}) « ri-^|A({z})|2j (1 - T - \zi\\..\zM?) . (85) 

provided < \zi\ < 1 for all eigenvalues, and zero otherwise. Here < T < 1 is the 
only non-zero eigenvalue of the deviation matrix. 

Eq. H85|) can be used to extract all n-point correlation functions defined as in Eq. ( 
lllf) . To achieve this it is convenient to use the Mellin transform with respect to the 
variable ( = 1 — T: 

POO 

Rnis;{zU)= dCC-'ia-O^'-'Rniizjn)] ■ (86) 

Jo 
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It is easy to notice that such a transform brings 'P{{z}) to the form suitable 
for exploitation of the orthogonal polynomial method The corresponding 

polynomials Pk{z) = z^\Jk + s are orthonormal with respect to the weight /(z) = 
|_j,|2(s-i)^^ inside the unit circle \z\ < 1. Following the standard route we find: 



= 



s{s + l)...{s + N -I) 



(87) 



where the kernel is 

K{Z1,Z;) = (/(Zi)/(Z2))^/' 



N-1 



S(t){x) + X — <j){x) j \x=ziz-i 

and (j){x) = {x^ — l)/(a; — I). Thus, the expression Ea. (|87(l can be rewritten as: 



Rn{s; {z}n) OC 



\2{s-l) 



s{s + l)...{s + N -1) 



-^s'g,({z}„); 



1=0 



qo{{z}n) = det 



lij' — l,...,n 



qn{{z}n) = det |ij = l,...,n 

and can easily be Mellin-inverted yielding finally the original correlation functions in 
the following form: 



Rni{z}n)\\.„\<l (X rl-^0(T -l + a)J2 llii^h 



1=0 



da 



1 



1 - T 



N-l 



(89) 



where 6{x) = 1 for a: > and zero otherwise. This equation is exact for arbitrary N. 
Let us now investigate the limit N ^ n>l and use, to the leading order: 



-a 



da 



1 



a 



N-l 



a~^ 



1 



a 



N-l 



This allows one to rewrite the correlation function as a determinant: 



Rni{z}n) OC T 



1~N 



a 



1 



1-T 



a 



e{T-l + a) 



(90) 



/N{l-T) , ^^ , 

i,j=i,...,n \a — 1 + 1 dx J ^ 

Further simplifications occur after taking into account that the eigenvalues Zi are 
expected to concentrate typically at distances of order of 1/N from the unit circle. 
Introduce new variables yi,6i according to Zi = {l — yi/N)e^^^ and consider yi to be of 
the order of unity when N oo. As to the phases 9i their typical separation scales as: 
di — 9j = 0(1/N). Now it is straightforward to perform the limit N oo explicitly 
and bring Ea.f I9()|l to the final form: 



i?„({z}„) oce-s det 



rfA(A + .g)e-5A(^(«.-fl.)-»(y.+w)) 



ij — l,...,n 
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with g = 2/T- 1 + 

The expression above coincides in every detail with that for rank-one deviations 
from Hermiticity, see Eas. (|42l43|l provided one fixes before the rescaled coupHng 
constant g and then remembers that the mean hnear density of phases Oi along the 
unit circle is v = l/(27r). The non- unitary matrices provide therefore an alternative 
way to address analytically " local- in-spectrum" universal features arising in the theory 
of chaotic scattering. Sometimes to deal with subunitary matrices is technically more 
advantageous than with their non-Hermitian counterparts. As an illustration of this 
statement we use subunitary matrices for extracting the distribution /(F) of the width 
r = 2Imz of the most narrow resonance among n — W/ IS. >> 1 falling in a narrow 
window [E — W/2,E + W/2] in the vicinity of a given point E in the spectrum. 
The latter quantity is of great interest in the theory of random lasing since it is 
related to fluctuations of the lasing threshold. The functional form of the distribution 
was found in those papers by employing plausible qualitative arguments of virtual 
statistical independence of the widths of the neighbouring resonances. For the simplest 
single-channel case we are able to show analytically that |57| 

The general case can be treated very similar. The distribution Ea.f l91(l is exactly that 
suggested in [SJ, but with renormalized effective coupling g = 5^(7 + 7"^) replacing 
the combination ^jl'Kv. We see that the two expressions coincide only within the 
weak coupling limit 7 — > 0, whereas the difference amounts to the factor 2 in the 
exponent for the case of perfect coupling 7 = 1. Let us briefly comment on the way of 
deriving the distribution Ea.f l91(l . Instead of extracting such a quantity from the joint 
probability density Eq. ( 131(1 we find it technically easier to consider its counterpart 
Ea.f I85|l and interprete the parameter N in Eq.JHHJ as the number n of resonances 
in the window. We know that in the limit n » 1 the eigenvalues Zfc = r^e^'' are 
situated in a narrow vicinity of the unit circle and their statistics is indistinguishable 
from that of the complex eigenenergies £, when the latter are considered locally, i.e. 
on the distances comparable with the mean spacing A. In particular, the distances 
1 — rfc from the unit circle should be interpreted as the widths of the resonances. To 
calculate the distribution Eq. ( I91|l we first notice that the form of Eq. ( I85|) allows one 
to integrate out the phases 9i by noticing that: 

/ fi.--/ ''^\{\r,e'''-r,e^''?^Ver{rl...,rl) (92) 
Jo 2^ Jo Svr 

where we denoted Per(a;i, ...,a;„) = ^1^ ■ ■ ■^n"^ ^'^d the summation goes over 

all possible permutations {a} = (ai,...,Q;„) of the set l,...,n (in fact here we deal 
with an object known as "permanent", hence the notation). In this way we come to 
the joint probability density of the radial coordinates only. Such a density written in 
terms of the variables Ri — rf has the following form 

Vt{Ri, ■■■,Rn) oc ri-^Per(i?i, ...,Rn)6 (1 - T - i?i . . . i?„) . (93) 

Identifying resonance widths with the distances 1 — r.^ , the probability distribution of 
the resonance ri closest to the unit circle: ri > r2,r3, ...,r„ is naturally described in 

+ The coupling constants Qa in this section are rescaled and correspond to ga/TTu(X) of subsection 
There is no difference aX X = 0. The coupling constants are always > 1. 
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terms of the function: 

Pt{Ri)= di?2.../ dRnVT{Rl,-,Rn) ■ (94) 

JQ Jo 

Again it is convenient to multiply this expression with T^^^ and to perform the McUin 
transform with respect to 1 — T. Then the integration over i?,— variables is simple to 
perform and after inverting the Mellin transform one finds: 



PT{Ri)^e{R^-l + T) 



1 



d 



n/2 



1 



1-T 
R^l 



(95) 



T"-! dRi 

The distribution Ea. l|91() follows from Ea. H95|l after rescaling i?i — l~2y /n, y = ttT / A, 
performing the limiting procedure n —^ oo and remembering A = 27r/n. 



3.4- Rank M > 1 deviations from unitarity 

Coming back to the general case we find it convenient to change: T^/^QiT^^^ 



Qi 



and separate integration over eigenvalues and eigenvectors of matrices Qi. The 
latter can be performed in a recursive way I — > / + 1, with the multiple use 
of the Itzykson-Zuber-Harish-Chandra Ea. (|28l) formula. After quite an elaborate 
manipulation, demonstrated in detail in appendix A2 one finally arrives at the 
following representation: 

^((4) o. ... i-r"'^ n (^ - ^) (»o) 



det(l-T)nc,<c. iT,.-T. 



dXe'^ 



iTrf A 



Ci<C2 

N 



dTr. 



k=l 



where we defined the diagonal matrices of size M as: f = diag(Ti,. 
diag(Ai, Xm) and used the notations: Tc = In (1 — Tc) and 



/(a, A) 



M 



1 Ils{^q)i^q ~ As) 



, tm) , A = 



(97) 



The distribution Ea. (|96|l is written for 



\Zk\ 



< 1 for any k = 1 , . . . , and 



vanishes otherwise. It shows great similarity to the pair of equations (|25I30|) describing 
eigenvalues of matrices deviating from Hermiticity. The remarkable feature of the 
distribution Eq. (|96|l is that it allows for calculation of all n— point correlation functions 
for arbitary N,n,M with help of the method of orthogonal polynomials. Again, the 
particular case M — 1 is quite instructive and can be recommended to follow first 
for understanding of the general formulae outlined below. 
To this end, we write 



N 



N 



iA({z})p n /(i^ i^fci'' ^) ^ n ^'^(^) 



fe=i 



N 

E 



.^1 ^"(A) 



A) 



ij = l,...N 



where the constants A^„(A) are provided by the orthonormality condition: 
d2zz™-i(z*)"-i/(ln \z\\ A) = <5,„,„Ar„(A) , 



(98) 



(99) 



|z|2<l 



which yields after a simple calculation A^„(A) — TrHcli 
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Now, by applying the standard machinery of orthogonal polynomials |49| one can 
find the correlation function: 



i?„(zi, ... 
as given by: 

Rn{zi, ■■■,Zn) (xV dXe' 



iTrfA 



{N -ny. 



N 



[]7Vfc(A)det K{z„z*;\) 



k=l 



where the kernel K is defined as: 



1 ^ 

K{zi,Z2;\) = - y^det (zA + nl)(ziZ2)' 

n=l 



(100) 



(101) 



(102) 



and the differential operator V is just the expression in front of the A— integral 
in Ea. l|96|l . It is worth noting that we used a version of the kernel which is not 
symmetric with respect to its arguments zi,Z2. This choice is legitimate, though 
unconventional, and is dictated mainly by our wish to avoid square roots which may 
induce complications when working with complex-valued expressions. 

In principle, all A— integrations in Ea. ljlOHl can be performed explicitly and the 
resulting formulae provide the desired general solution of the problem. However, for 
arbitary N, M, n the results obtained in that way are still quite cumbersome. Some 
examples are given in ^E^' . In the theory of quantum chaotic scattering we, however, 
expect a kind of universality in the semiclassical limit. Translated to the random 
matrix language such a limit corresponds to TV — > oo at fixed n, M. Still, extracting 
the asymptotic behaviour of the correlation function i?„(zi, z„) from Ea. (|lUll) in 
such a limit is not a completely straightforward task. A useful trick is to notice that 
Ea. H101|l can be rewritten as: 

M 

.,Z„)CX Fq- 
qi = l,...,qn = l 

N 



Rk{zi, 



qi,...,q„ {{Tc\ z}) 



1...M 



{{T,;z})^B{T,}det 
d 



E 

,fc=i 



Cl<C2 



OTc, 



M 



n 



-d 
exp —iXcXc 



( -d 



\dTM 



+ k [ziZ*) 



fe-i 



(103) 
(104) 



K) 



where we used the notation: 
B{T,} = 



1 



M 



TT- 

nc,<c. (r,, -r,,)li (i-r,) 

Introducing now the auxiliary differential operator = Ilj^i 11^(^4^) (sr" ~ WT^ 

and considering its action upon the ratio -Fl3i,...,q„/iB{Tc} one can satisfy oneself that 
in the limit TV ^ M, n the leading contribution to i^5i,...,g„ is given by: 



M~N 



91, 



,9n 



M 
c=l 



N-AI 



n 1...M 



n n i^-w 



det 



if(z„z*;{fj) 
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where the kernel is given by 



N M 



k=l c=l 
M 



{N - M) 



fe-i 



(106) 



n 



Tr dx 



1 - X '^=^'^; 



Sq,,c\ri\zj\^ 



and we used the notation: Tc — 1 — exp ^Tc — X]J=i ' 

Further simphfications occur after taking into account that the eigenvalues Zi are, 
in fact, concentrated typicaUy at distances of order of 1/A^ from the unit circle. Then 
it is natural to introduce new variables yi,4>i according to Zi = (1 — yi/N)e'^^' and 
consider y.i to be of the order of unity when N oo. First of all, one immediately 
finds that: 

ji-^nf^l =-pLx:y.^i (107) 



c=l \ 



As to the phases 9i, we expect their typical separation scaling as: 9i — 9j = 0{\/N). 
Now it is straightforward to perform explicitly the limit iV — > cx) in Ea. (|106() . 
Combining all factors together, one brings the correlation function Ea. H103|l to the final 
form coinciding in every detail with that obtained for random GUE matrices deformed 
by a finite rank anti-Hermitian perturbation, see Eas. (|42l43|) . with gc = 2/Tc—l. One 
should only remember that the mean density of phases di along the unit circle is 
ly = l/(27r) and take this factor into account replacing the (semicircular) density 
factor v^x) whenever necessary after fixing the rescaled coupling constants ga- This 
completes the proof of equivalence, in the large N limit, of spectral properties for 
finite-rank deviations from Hermiticity and from unitarity. 

We have seen already that the correlation function of moments of characteristic 
polynomials played a very important role in the theory of non-Hermitian matrices. It 
is natural to consider similar objects for general subunitary matrices A as well. So far 
only the simplest objects of this kind were calculated Wl : 



Iizi,Z2) 



det (zil - a) det ( zXl - 



(108) 



where the angular brackets stand for the averaging over the probability density in 
Eq.Q. We present the final result, which has the form of a Hankel determinant, in 
terms of eigenvalues Gi = 1 — of G: 



I{ZUZ2)^1[ 



{N 



j!(j + l)!(7V + j)! 



det [f, 



(109) 



N 

E( 

fc=0 



{N + i + j - fc)!(2n -2 + k-i- j)! 
{N + 2n- 1)! 



l<il< 



E 

2<...<ik<N 



■G, 



The particular case n = 1 of the above formula was presented in In fact, the 

expression Ea. (|l()9(l is valid for any Gi > 0, being not at all restricted to the subunitary 
case < Gi < 1. For random unitary matrices all G^ = 1 and one can show that the 
expression above are indeed in agreement with the known result for the moments of 
characteristic polynomilas of unitary matrices |82[ . 
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4. Statistics of Non-Orthogonal Eigenvectors 

Let us denote \Rk) and {Lk\ to be the right and the left eigenvectors of the matrix H 
corresponding to the eigenvalue Zk- This means 

n\Rk)=Zk\Rk) , {Lk\n^{Lk\zk (110) 

\Lk) = zl \U) , rC = (Rk\ zl . (Ill) 

These eigenvectors form a complete, bi-orthogonal set and can be normalized to satisfy 

(L„i I i?„) = 6mn (112) 

The most natural way to characterize the non-orthogonality of the eigenvectors is to 
consider the statistics of the overlap matrix — {L,-a \ L„) {Rm \ Rn)- This matrix 
naturally enters many calculations that operate with non-Hermitian Hamiltonians, 
e.g. in the description of the particle escape from the scattering region ("norm 
leakage" j84|'). Moreover, the entries 0„m have a direct physical meaning for lasing 
media: diagonal ones are the so-called Petermann (or " excess noise" ) factors accessible 
experimentally, whereas the off-diagonal entries represent cross-correlations between 
the thermal or quantum noise emitted into different eigenmodes |6(J) . 
Thus, two correlation functions 53 : the diagonal one 

0{z)^I^Y.^"-Siz~Zk)\ (113) 

and the off- diagonal one 

0{z, ^"-^^ - ^ - ) (114) 

are natural characteristics of average non-orthogonality of eigenvectors corresponding 
to the resonances whose positions in the complex plane are close to the complex energy 
z. Here 5{z) stands for the two-dimensional (5-functions of the complex variable z. Let 
us note that for any ensemble with orthogonal eigenvectors and complex eigenvalues 
z (e.g. for the so-called normal matrices) one obviously has 0{z) equal to the mean 
density of complex eigenvalues. For the same situation the off-diagonal correlator 
vanishes: 0{zi,Z2) = 0. 

Both diagonal and off-diagonal eigenvector correlators were introduced and 
explicitly calculated for the strongly non-Hermitian (Ginibre) ensemble of non- 
Hermitian matrices by Chalker and Mehlig For the ensemble J'n = H ~ iT 
pertinent to chaotic scattering both types of eigenvector correlators were found for the 
regime of very strongly overlapping resonances when widths typically much exceed the 
mean separation (^HI ■ We already discussed that such a regime corresponds to a large 
number M ^ 1 of open channels, and a kind of self-consistent Born approximation 
(or equivalent approximation schemes, see 55 ) is justified in this case. A general 
non-perturbative expression for the diagonal correlator 0{z) can be extracted from 
the following relation |55j : 

0(z) = -lim( f Tr . | ) . (115) 

Although a direct calculation of the right-hand side for weakly non-Hermitian matrices 
is difficult, Frahm et al. ,54, managed to get to the result valid for any number 
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M of open channels by employing a kind of heuristic analytic continuation scheme 
suggested in JT] and known to reproduce the exact expression for the resonance 
widths distribution, like that in Ea. H74|) . Therefore we have good reasons to expect 
the result for the eigenvector correlator found in 54' 

W^O(.)M.) = l + -^-^M_ , ny)-~jy^TA^)iTA:c) (116) 
where y — and 

should be valid as well, though it is still desirable to verify this formula independently. 
No non-perturbative results for the off-diagonal eigenvalue correlator 0{zi,Z2) were 
reported so far, to the best of our knowledge. 

In a recent paper |57j exact non-perturbative expressions for both diagonal 
and off-diagonal eigenvector correlators were provided for the case of a system with 
broken time-reversal invariance coupled to continuum via a single open channel. The 
expression for the diagonal correlator was found to agree with the formula Eq. H116|l and 
was verified by comparing the analytical expression with the result of direct numerical 
diagonalization of the matrices JTat. In fact, numerically it is easier to compute 
smoothed averages, such as the mean number of eigenvalues {n{Lx, Ly)) inside the 
rectangle A — —Lx/2 < Rez < Lx/2 , < Imz < Ly ^. Such a quantity can 
be obtained from the mean density p{z) by a simple integration over the rectangular 
domain A. Similarly one can define the function Oi{Lx, Ly) as the integral of the 
diagonal correlator 0{z) over the same domain. Numerically the latter quantity should 
be compared with J^z^eA ^kk- 

The formula for the off-diagonal correlator 0(zi, Z2) found in 57 is given by: 

0{z,,Z2)^N(^^y e-a(y^+y^) I' dX, C rfA2(.g + Ai)(.g + A2)e^"(^^+^^) (118) 



Xe-!'2(^i-^2)/2 



-1 J-i 



where we introduced Imzi 2 = ^1,2 = ^ -\- assumed that 57 ~ A and denoted 
oj = TrfJ/A. Again, for numerical reasons it is easier to calculate the smoothed 
average ^ Omn which should be compared with the corresponding integral of 
the function 0{zi,Z2)- 

Let us now outline the derivation of the formulae Eas. (|116lll8|l . The main idea is 
to use the fact that the complex eigenvalues Zk of Tie/ / = J^n are resonances, i.e. the 
poles of the M x M scattering matrix S{£) in the complex energy plane £. Denoting 
by y a (nonunitary!) matrix of (right) eigenvectors of Ti-eff = VZV^^ we have, in 
particular, the identities: 

K \ ■' / J nni ^ J nm 

Let us also note that according to the definition the entries of the overlap matrix Omn 
are given in terms of V as 0„„, = {^^^}„„ [v^'^ [^"^]^} • 
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121) 



Substituting the relation E — Ti-efj — V{E — Z)V~^ into the formula for the 
scattering matrix Eq.(|3Jl we can extract the behaviour of the scattering matrix 
elements when the (complex) energy parameter £ approaches an individual pole: 

Sab ^ z„) « -2^^ (W^V)^^^ {V~'WU-^ (120) 

A similar formula holds for the Hermitian conjugate of the scattering matrix: 

These two equations immediately yield: 

Tr{5(£i^z„)[5t](£2^z„0}= (122) 
(£i-z„)(^2 -<) ^ ^ ^ ^""^ 

Now we notice that 2TTiWW^ = Hef f —Ti-lj: j and further exploit the relations Ea. (|119() . 
This immediately yields the following relation between the trace of the S-matrix 
residues and the overlap matrices of the left and right eigenvectors: 

Tr [Res5(f )Res S^S*)] - « - ^n) (z„ - z*J 0,nn ■ (123) 

Such a relation is of general validity, but for an arbitary number of open channels 
M it seems to be of no obvious utility, due to difficulties in evaluating the ensemble 
average of the trace of the residues on the left-hand side. However for the particular 
case of a single open channel M — 1 the scattering matrix coincides with its 
determinant and therefore has a different representation as (see Eq.Q) 

S{£) = 1[t^^ hence 5t(f*) = ]J (124) 

k=l k=l 

In fact, this formula is model-independent and follows, up to an irrelevant "non- 
resonant" phase factor, from the requirement of S'-matrix analyticity in the upper 
half-plane and unitarity for real energies. This expression substituted in Ea. H123|l 
yields immediately the relation: 

(z —Z* Y Z„ ~ Zh Z* ~ Z*i 

giving the eigenvector non-orthogonality overlap matrix in terms of the complex 
eigenvalues Zk- This provides us with a possibility to find the diagonal and ofi^-diagonal 
correlators, Eas. (|113lll4|l . by averaging Omn over the known joint probability density 
of complex eigenvalues for the single-channel scattering system Ea. (pTT|l . Indeed, one 
may notice that 

- N-2 

0{z) = ^e-M-^7^-(-^-i)T.]e-f (^^+^-) (det[(z - Jt) _ ^)])^^_^ (126) 

where J^n-i stands for the non-Hermitian matrix of the same type as jTw but of the 
lesser size {N — 1) x {N — 1), and with coupling 7 replaced by a modified coupling 
7i = 7 — Imz. Analogously 



X <det[(zi - J^) (z* - J^) {Z2 ~ J) {zl ~ J)]) 



J'n - 
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where Jn~2 is of the size (A^ — 2) x (A^— 2), and with coupUng 7 replaced by a modified 
coupUng 72 = 7 — Imzi — IniZ2. 

In this way the problem is again reduced to calculating a correlation function 
of the characteristic polynomials of large non-Hermitian matrices as that in Eq. (|33|) . 
The scaling limit iV > 1 such that Im2:i,2 = ri^2 ^2^^ Re [zi - 22) A cx N^^ of 
the resulting expressions just yields the formulae Eas. (|116lll8(l above. 

5. Conclusions and perspectives 

In the present paper we discussed properties of random matrix ensembles 
which can be naturally called non-Hermitian/non-unitary deformations of classical 
Hermitian/unitary ensembles. 

Guided by the example of resonances in quantum chaotic scattering, we revealed 
the existence of a non-trivial regime of weak non-Hermiticity. The regime can be 
defined as that for which the imaginary part IvciZ of a typical complex eigenvalue 
is of the same order as the mean eigenvalue separation A for the corresponding 
Hermitian counterpart H. We systematically demonstrate that for such a regime 
the correlation functions of complex eigenvalues describe a non-trivial crossover from 
the Wigner-Dyson statistics of real eigenvalues typical for Hermitian/unitary random 
matrices to the Ginibre statistics in the complex plane typical for ensembles with 
strong non-Hermiticity. The latter is defined by the condition < Trff^ >(x< TrF^ > 
when N 00. Another important new feature of the non-Hermitian matrices is the 
emerging non-orthogonality of eigenvectors, and we discussed a few results available 
in the literature. 

Among challenging problems for future research we would like to mention 
extending our understanding to other symmetry classes of weakly non-Hermitian 
matrices, most importantly to complex symmetric and real asymmetric cases. In fact, 
the joint probability densities of complex eigenvalues for some particular cases are 
known, see e.g.^HZHI, but no progress in extracting any correlation function beyond 
the mean eigenvalue density was reported so far. For non-perturbative calculations 
(replica trick) of the mean eigenvalue density for various non-Hermitian ensembles 
see the recent paper |74) . A consistent non-perturbative treatment of non-Hermitian 
matrices with "chiral" symmetry is still an open problem as well. Next to nothing 
is known about statistics of non-orthogonal eigenvectors for other classes of non- 
Hermitian matrices, in particular, for the important case of the time-reversal invariant 
scattering. 

Our considerations were mainly restricted to a Gaussian probability measure for 
the Hermitian counterparts of the considered matrices. An interesting open issue is 
the universality of the obtained results versus modification of the probability measures 
for Hermitian and anti-Hermitian parts. Relatively little is known in this direction 
and the available results are restricted by some specific types of the probability 
measure |44) . Within the non-linear cr— model formalism one can demonstrate a 
kind of universality by non-rigorous, heuristic methods for matrices with independent 
entries. A challenging problem is to elevate our understanding of properties of almost- 
Hermitian random matrices to the level typical for their Hermitian counterparts. 

Some interesting questions still remain to be answered even for non-Hermitian 
deformations of GUE matrices. To ensure their applications to the theory of chaotic 
scattering, we first have to improve our knowledge on eigenvector statistics. In 



Random matrices close to Hermitian or unitary 



38 



particular, one needs to study the case of more than one channel and the problem 
of understanding fluctuations of the non-orthogonality overlap matrix Omn- 

From all these points of view a most perspective direction seems to be the method 
of reducing various correlation functions of interest to products (and, in the general 
case, ratios!) of spectral determinants (characteristic polynomials). For Hermitian 
random matrices closely related objects were recently quite intensively treated by 
various analytical techniques |8(ij and the issue of universality proven to be amenable 
to a rigorous mathematical treatment [HZI ■ 

Recently, there was considerable interest in understanding interplay between the 
weak non-Hermiticity and Anderson localization 76 and in statistics of the resonances 
in systems with localized eigenfunctions Clearly, this issue deserves further 

attention. 

All questions about non-Hcrmitian matrices can equally be asked for matrices 
deviating from unitarity. Actually, some statistical properties of general subunitary 
matrices were under investigation recently as a model of scattering matrix for systems 
with absorption, see jHH- One may hope that this class of random matrices will find 
growing applications in the future. 

6. Acknowledgements 

The authors are much obliged to B.A.Khoruzhenko, B.Mehlig, M.Titov and K. 
Zyczkowski for collaboration on different aspects of non-Hermitian random matrices 
discussed in the present paper. Communications and discussions with many 
colleagues over the years, especially with G.Akemann, Y.Alhassid, A.Garcia-Garcia, 
S.Grebenschikov, T.Guhr, I.Guarneri, F.Haake, G.Hackenbroich, T.Kottos, J.Keating, 
V.Mandelshtam, L.Pastur, D. Savin, H.Schomerus, T.Seligman, V.Sokolov, T.Wettig, 
H.-A.Weidenmueller, and J. Verbaarschot on various questions related to the subject 
are much appreciated. 

This work was supported by SFB 237 " Unordnung und grosse Fluktuationen" and 
by EPSRC grant GR/13838/01 "Random matrices close to unitary or Hermitian" . 

Appendix A. Various integrals over the unitary group 

Appendix A.l. Itzykson-Zuber integral for the case of matrices having reduced rank 

Denote Uc , c = 1, iV columns of the unitary N x N matrix U G U{N). They are 
orthonormal A'^— component complex vectors: uJ^Ucj — SciC2- To evaluate the integral 
H28|l for the case of a matrix F of reduced rank M < N (i.e. when only M eigenvalues 
denoted by 71, ... , jm are nonzero, the rest N — M being zero) we notice that the 
combination Tr K WTU — J2c=i uJiCuc obviously depends only on the vectors Uc 
with 1 < c < M. These column vectors form a. N x M matrix which we denote Vm, 
and the (Af— component) rows of this matrix we denote as v|| , j = 1, ...,N. Then 



to the condition of the matrix /a/(v) to be just the identity matrix: /Af(v) = 1^/. 
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Therefore the integration measure over the manifold of matrices Vm must be given by 

N / N \ 

d \Vm] oc W d^w, 5\Y^^r,® v] - 1m (A.l) 



j=i v=i 

N 



dF, 



M 



„*Tr#M(i:f=iV,«5vt-i„) 



where we used the (matrix) Fourier-representation for the (matrix) S— function, with 
Fm being the manifold of Hermitian M x M matrices. 

These facts allow us to rewrite the IZHC integral for the present case as: 

N 



r ^ -1 

cx / d/'Afe-^TrfM -Q (^Fm - ^m) ■ 

These integrals to be well-defined we assume a positive infinitesimal positive imaginary 
part added to Fm in the denominator. Now we introduce the Hermitian matrix 
F by Fj^i = jj^j F^j^j which yields the Jacobian factor: dFM = [dei^M] dF. 
Diagonalizing F — U^/AmU m we use the matrix of eigenvalues Km = diag(Ai, Am) 
and M X M unitary matrices U m as new variables of integration, with the usual change 
in the measure as dF = I^^ {kM)dkM[dU]M ■ After these manipulations we arrive at 
the relation 

r r N M 

/ [dUU e-TrKC/trc/ ^ ^^^.^^^^f-^ / dAMA2(AM) W W — ^ / e-Trw^AM^/M 

c — 1 j — 1 

where we used the standard IZHC formula Eg. pSfl for the integration over Um- 
Finally, at the last step we use the symmetry of the integrand with respect to 
permutations of the arguments Ai, Am and arrive at the required extension of the 
IZHC integral (P|) : 

/J J- M—N f ^ 
[dU]i, e-^Tri^r/tm ^ j dkM^M) n /c(^c), (A.4) 



where 



\Aj=i (Ac - kj) 

which is valid for the case when one of the matrices is of a reduced rank M < N. Let 
us finally satisfy oneself that when both matrices K and F are of the full rank: M = N 
our expression Ea. l|A.4|) reproduces exactly the familiar Itzykson-Zuber formula 128() . 
To this end we use the Cauchy identity: 



det 



1 \ _ A(A)A(if) 
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which aUows us to write the write-hand side of (|A.4() as 
AIWT / rfAe- S£i 7,A, det (^) 



cj=l,Ar (A.5) 



where the suramation goes over aU permutations T — (kj-^, kj^^) of the set 
(fei, fcAf), with the factors (—1)"^ standing for the signs of permutation. Each 
integration yields a factor: iex'p{—i'yckj^). Taking into account the factors (— 1)"'" 
allows one to present the result as a conventional Itzykson-Zuber determinant. 

Appendix A. 2. Integration over unitary group for subunitary matrices 

We start from Ea. l|83|) and want to perform all integrations over matrices Qi and 

complex vectors a. Rescaling all variables according to Qi = \fi' Qi with 

Qo = and Qn — and a = a and taking into account all Jacobians from 
the transformations we arrive at 

V{{Z))^\l^{{Z})\^ YldQ^H d^kS{Q,-Q,^i-a(g>a.^) 
i=i i=i"' 

N 

[]<5(|zi|2...|z,p-det[l-Q,]) detf^^-^. (A.6) 



N 

X 

1=1 



The following strategy is to diagonalize all matrices Qi as Qi — UiXiUi ^ with 

j<k I 

and then perform the integrations over the eigenvectors ?7i. We start with j d^a.6{Qi — 
a (g) a'f ) . Since we integrate over all directions of the complex vector a the integral 
does not depend on Ui, i.e. the eigenvectors of Qi and is given by (cf. Ea. l|A.l|) : 

J d^a (5(gi - a(g) a^) cx J dfl expiiTrflQi) / dct{e + ifl) 



Here we diagonalized the matrix f2 and integrated out the corresponding eigenvectors 
with the use the IZHC formula Ea. l|28|) . with e > (cx unity matrix) being an 
infinitesimal regularisation matrix. The result is a complicated expression in terms of 
^-functions and step functions. We now proceed to calculate the next factor in HA.6() . 
Using again the IZHC integral two times we obtain 

d,m J d-a SiQ, - - a « at) cx ^J^U^^^^^^^^^ (A'^) 

We used here that (|A.7|) was independent of [/i and we see that (|A.8|I is independent 
of U2- Using the latter fact on the next stage, we are able in this way to perform the 
integrations over all Qi-diagonalizing matrices Ui. There remain the integrations over 
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all the eigenvalues A^-"''' of the matrices Qi. The resulting expression for V{{Z}) may 
now be written as 

V{{Z}) cx ^ ___| A({Z})|2 j W dQ, 0{f - QN-i)...0iQ2 - Ql)[Al^ 



N 

l[5{\z,\\..\zi\'-det[l-Qi]) . (A.9) 



N 

X 
1=1 

with diagonal M x Af matrices Qi, and the notations o?Qi := YifLi '^■^i^ 

M M 

^(02-Qi):-n^(^2''--^i'^) ' det(l-Qi) :=n(l-Ai*^) (A.IO) 

i=l 1=1 
l...Af / Q /9 \ *^ 

An important simplification occurs if we go to logarithmic variables = ln(l — Qi) 
(again these are diagonal M x M matrices with > n > ln(l — T) = rjv )• We have 

5 - det [1 - Qi]) 9{Qi - g,_i) = M E 1^ I^jI' - E 1^(1 - 



det(l - QO 

A/ 

X Y[ 9{\n{l - A|i\) - ln(l - A^^)) 



j=i k=i 



and there is the special relation 

1...A I / ^ „ \ Af 

j<k V^^l / (=1 

Taking into account the Jacobian from the transformation our joint density takes the 
form 

" Tm T / n '^-^ 9{n - r,)...e{r^^, r^) 

det(l - T) {{^^kiTj - Tk) J 

N 

l[5{\n\z,f + ...+ln\z,f-TrTi) . (A.12) 



\.j<k\-^3 

N 

X 
(=1 

Now we see that the operator Ai can be shifted by partial integration, it commutes 
with Trri and then acts as A2 on 0{ti — T2). This argument can be repeated, such 
that finally the operator A^r appears in front of the remaining integral. Changing 
again the integration variables we obtain 

V{{Z}) cx dctr^^-^|A({Z})p ^ f ^ . e{-Ti)5{\n\zi\^ -Ttti)] diY t,^t) (A.13) 
det(l-T)n,<,(T,-T,) yfJL I K V K \ i\ i)\ ^Z. ' M 

with T := ln(l — T) and A :— nj <fc^ ( a^TT ~ drW ) ' Writing down the Fourier 
representations of the ^-functions in (|A.13I) the integrations over the diagonal M x M 
matrices can mainly be performed yielding 

/ l[[dnS{\n\zi\^-TTTi)]s{j2n-T] = ^ expHTrcDr)n/(lnki|','^)(A.14) 

J{r^<0}l^l \^=l I -I L^'^J , = 1 
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with 

^ / dneM^ ^ y exp(^a;,a) 

and with an infinitesimal e > 0. The integration over can be performed in the 
complex plane yielding the restriction a < 0, i.e. in (|A.13p all Zi lying inside the 
unit circle (In < 0). Thus we have derived the full joint density of eigenvalues of 
subimitary matrices (|95|l . which can be used to derive all correlation functions. 



Appendix B. Derivation of Ea.()32[l 

Notice that in the limit N ^ oo we can replace N in the exponent in Eq. (|24() by 
N — n. Then, similar to the Gaussian case, see Ea. H22|) . we first write: 



X n ni^. -^'l' /rfAa(A,7)A(A) I] E -fT 
1=1 i=i 1=1 c=i f^}^ 



(Ac — As 



where 



a(A,7) = e-'i:^".-^=nE-ff7r-AT= ^ ■••E ^'^(^,7) (B.2) 

1 = 1 K = l K„ = l 

p-ii:fll7cAc 

^k(A,7) 



n M 

n n (A., -A, 

] = 1 s{^Kj) 



In Eq. (|B.2(I we denoted k — (ki,...,k„) and 7c = 7c — '^c.kj Imzj, where 

Sc,K, — I a c ~ Kj and zero otherwise. This allows us to write: 

RniM)^'- ^— ;f^^E---E^«(7,W), (B.3) 

K.i — 1 Kji — l 

where 



F.(7,{z}) = i?(7)/d'{ae" ' ' I AlUDPH 11 1^^" " I' (^.4) 

1=1 j=i 

rfA A4A,7)A(A) n E-ff 

i = l c=l j-}^ , 



(Ac - As) 
and we denoted 



i?(7)= 



e 2— i-c=l 7c 



A(7)det7^-"-^^' 

Now we introduce the differential operators = YYj=i Y[t{^K ) ~ Wf 

and consider its action on the ratio Ff^{'-f) / B (j) from Ea. ljB.4|) . On one hand, we 
notice that: 
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and, as a result the action of Z)^ on Ff^{'y)/B{j) brings the ^-integrand exactly to the 
form that allows us to rewrite it as an average of a product of determinants (cf.Eq.|221)). 
The averaging goes over the density of joint probability of eigenvalues ^i, ■^2, ■■■,^N-n 
of a non-Hermitian matrix J7v-n(7) = Hjsi-n + i^ , with nonzero eigenvalues of F being 
71, ■■•,7m: 



A(7) dety 



J]^ det [zj - J]s-n{l)] 

\3 = 1 



(B.5) 



GUE 



On the other hand, to ensure existence of a well-defined limit ^ oo of the 
correlation functions Rn{{z}) in Ea. ljB.3|) the quantities F^{'-^) must in that limit 
behave regularly as functions of 7c, i.e. F„(7)|Ar^i ^ CNFj^'\^) where Cn does not 
depend on 7 and ^^"■'(7) has a finite limit when iV ^ 00. Because of this property in 
the limit N ^ n the action of on 



F,(7)/S(7) = F,(7)A(7)det7^-"-*'£ 
to the leading order amounts to: 



[F,{j)/B{^)]^F^{j)A{j)D, 



det 7 



N-n~M 



e 2 



(B.6) 



Indeed, each 7-differentiation of the terms in the square brackets in Ea. l|B.6|) brings a 
factor of the order of N, which is much larger compared to the result of differentiating 
the factor i^„(7)A(7). Performing the remaining differentiations explicitly and 
considering N — n^N — n — MRiN, we find: 

[F47)/S(7)] \N»n,M « (27V)"(*^-i)^^«(7)A(7) 

n M 



X det 7 



e 2 



n n - 9s 

where we introduced the notation gc = ^(7c + 1c 
Comparing this expression with Ea. (|B.5p yields 

. -I N~n~M 



J^k(7,{2})« 



C^{{z}) 



A(7) 



n;=in.U)(5..-ff.)A(7) 

where we denoted: 



det 



N-n ^ M (2 ~2\ 

e — ^5:„i(7e-7c) (B.7) 



) ■ (B.8) 

' GUE 

In the limit when N 00 and M finite, almost all eigenvalues of J have imaginary 
part of the order <C 7c IHOl- Rescaling the imaginary parts yj — Nlmzj one finds 
that 



A(7) 
A(7) 



det 



7 



N-n-M 



in the limit n, M <ti N ^ 00. As will be clear later on, in the same limit one can 
substitute 7 for 7 in the GUE averaged product of determinants. Combining all the 
factors we arrive at Ea. H32|l . 
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